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ABSTRACT 



We have analyzed a nonsingular model with a variable cosmological term following the Carvalho et al. 
ansatz. The model was shown to approximate to the model of Freese et al. in one direction and to 
the Ozer-Taha in the other. We have then included the effect of viscosity in this cosmology, as this 
effect has not been considered before. The analysis showed that this viscous effect could be important 
with a present contribution to the cosmic pressure, at most, of order of that of radiation. The model 
puts a stronger upper bound on the baryonic matter than that required by the standard model. A 
variable gravitational and cosmological constant were then introduced in a scenario which conserves 
the energy and momentum in the presence of bulk viscosity. The result of the analysis reveals that 
various models could be viscous. A noteworthy result is that some nonsingular closed models evolve 
asymptotically into a singular viscous one. The considered models solve for many of the standard model 
problems. Though the introduction of bulk viscosity results in the creation of particles, this scenario 
conserves energy and momentum. As in the standard model the entropy remains constant. We have 
not explained the generation of bulk viscosity but some workers attributes this to neutrinos. Though 
the role of viscosity today is minute it could, nevertheless, have had an important contribution at early 
times. We have shown that these models encompass many of the old and recently proposed models, in 
particular, Brans-Dicke, Dirac, Freese et al, Berman, Abdel Rahman and Kalligas et al. models. Hence 
we claim that the introduction of bulk viscosity enriches the adopted cosmology. 
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Chapter 1 
Introduction 



The standard model was successful in describing the Universe, since it predicts the existence of the 
observed cosmic microwave background radiation having a temperature of 2.75K. The other success was 
the formation of light elements during the first few minutes after the big bang. However, this model is 
fraught with some vexing problems. These problems remain unsolved within the context of the standard 
model. An extension to this model becomes inevitable, this extension may include an extension to the 
theory of general relativity. 

Einstein introduced a cosmological term to the field equations to obtain a static Universe, but later, 
following the discovery of cosmic expansion, discarded this term. However, there are good reasons to 
believe that this term had an appreciable contribution at the beginning of the Universe. From a point of 
view of quantum field theory this term could correspond to a vacuum energy which was found to be very 
much bigger than the observed value. This discrepancy is known as the cosmological constant problem. 
In Particle Physics, the solution of this problem amounts to finding a mechanism that requires this term 
to vanish. An alternative to this view is to assume that as the Universe evolves, this term evolves and 
decreases to its present value. According to this view the cosmological term is small because the Universe 
is too old. Recently, cosmological models with this term decreasing with time were proposed (Ozer & 
Taha, Freese et al, Gasperini, Chen & Wu, Abdel Rahman, Carvalho et al. ). These models are known 
as vacuum decaying models. 

In Chapter 2 and 3, we review the standard model concentrating on its successes and shortcomings. 
Chapter 4 deals with a general review to the nonstandard cosmology. A nonsingular model depending 
on the Carvalho et al. ansatz (A = 3/3if 2 + where a, (3 are arbitrary dimensionless constants) is 
introduced in Chapter 5. In Chapter 6 we present a closed nonsingular viscous model with a cosmological 
constant of the Chen and Wu type and discuss its implications for the evolution of the Universe. Chapter 
7 is devoted to a singular model with varying cosmological and gravitational constants containing bulk 
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viscosity varying as p n , where p is the density of matter and n is a positive number. Various models 
are reproduced with appropriate values for n. A flat viscous model with (5 = is presented in Chapter 
8. The result of the analysis shows that G increases with time in both radiation dominated and matter 
dominated phases of cosmic evolution. These results are in accordance with the previously found results 
of chapter 7. An appendix followed by some published papers is presented at the end of the thesis. 

Finally we remark that in these cosmologies we have not discussed the issue of influence of bulk 
viscosity on the formation of galaxies. We will treat this in a future work. 
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Chapter 2 

The Standard Model 



2.1 Background 

A cosmological model is a model of our Universe, taking into account and using all known physical laws, 
predicts (approximately) correctly the observed properties of the Universe, and in particular explains in 
detail the phenomena in the early Universe. In a more restricted sense cosmological models are solutions 
of Einstein's field equations for a perfect fluid that reproduce the important features of our Universe. All 
cosmological models which differ only near the origin of the Universe must be accepted as equally valid. 
In fact a series of solutions are known which are initially inhomogeneous or anisotropic to a high degree, 
and which then increasingly come to approximate the observed Universe. All cosmological models which 
yield a red-shift(z) and a cosmic background radiation can hardly be refuted. The possibility cannot 
be excluded that our Universe is neither homogeneous nor isotropic, but has those properties only 
approximately in our neighborhood. Every model is of course a great simplification of reality, and only 
by the study of many solutions can one establish which simplifications are allowed and which assumptions 
are essential. 

2.2 Einstein's field equations 

These equations describe the gravitational field resulting from the distribution of matter in the Universe. 
They are nonlinear partial differential equations. 

?V - \ n 9»v = 8ttGT^ u (2.1) 

where R^ u is the Ricci tensor and 1Z is the curvature scalar, T^ v is the energy-momentum tensor of the 
source producing the gravitational field, and G is Newton's gravitational constant. 
The energy momentum tensor satisfies the following requirements. 
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(i) Ijjj, is symmetric with respect to interchange of ji and v . 

(ii) T^y is diver genceless, for energy and momentum to be conserved, 
(Bianchi identity): 

T> xv ;v = Q (2.2) 

The Ricci tensor is defined by 

_ d 2 ln^ dT%, _ e dln^ 

and 

n = gTRua, , (2.4) 

where g^ v is the metric tensor, g its determinant and T l - k are the Christoffel symbols related to g^ v by 

T )k = 9 ta (9aj,k + 9 a k,j ~ 9jk,a)- (2-5) 

2.3 Robert son- Walker Metric 

The distribution of matter and radiation in the observable Universe is homogeneous and isotropic. While 
this by no means guarantees that the entire Universe is smooth, it does imply that a region at least as 
large as our present observable Universe (Hubble's volume) is smooth. So long as the Universe is spatially 
homogeneous and isotropic on a scale as large as Hubble's volume, for purposes of our local Hubble's 
volume we may assume that the entire Universe is homogeneous and isotropic. This is known as the 
Cosmological Principle. The metric for the space with homogeneous and isotropic sections which being 
maximally symmetric is the Robertson - Walker (RW), 

ds 2 = dt 2 - R 2 {t){ 2 + r 2 (# 2 + sin 2 Od 2 ^)) (2.6) 

where (t,r,9,4>) are coordinates (referred to as comoving coordinates), R(t) is the cosmic scale factor, 
and with appropriate scaling of the coordinates, k can be chosen to be 1, -1, or for spaces of constant 
positive (closed), negative (open), or spatial (flat) curvature, respectively. The time coordinate in the 
above equation is the proper (or clock) time- measured by an observer at rest in the comoving frame, 
i.e. (r, 6, 4>) = constant 



2.4 Expansion of the Universe 



The RW line element contains the function R(t), which can be any function of time, t. This line 
element necessarily requires that the Universe cannot be static. Because of the presence of the term 



R(t), an element of the spatial distance dl changes with time, i.e. the distance between any pair of 
galaxies changes with time. Because of homogeneity, the cosmological fluid can not sustain any pressure 
- gradient, the concomitant non-gravitational forces are absent. But when the only force present is 
gravitation, a static Universe is evidently not possible, it must either contract under gravity or expand 
against gravity. However, the analysis of the observed data from distant galaxies shows that the Universe 
is expanding. The recession velocity of a galaxy turns out to be related to its distance (d) by a simple 
law v = Hd where the proportionality constant H is the Hubble constant. Note that the recession has 
no effect whatsoever on individual bodies since a homogeneous medium generates no gravitational field 
inside a spherical cavity. 

2.5 Distance measurements in RW metric 

The determination of distance in astronomy is mostly done using the concepts and ideas of a three- 
dimensional Euclidean space. We therefore want to describe briefly how the laws of light propagation 
in RW metrics influence this distance. One possible way of determining the distance of an object is to 
compare its absolute brightness L, which is defined as the total radiated energy per unit time and is 
regarded as known, with the apparent brightness I of the energy reaching the receiver per unit time per 
unit surface area. The Luminosity Distance Dl is defined by 



D L = ^L/4ttI (2.7) 

so that in Euclidean space luminosity distance and geometrical distance coincide. eq.(2.6) can be written 

as 

ds 2 = dt 2 - R 2 [d X 2 + f(x)(d0 2 + sinfl 2 # 2 )] (2.8) 

where 

^ VT— kr 2 ^ ^ 

In the RW space-time there exists a complicated relationship between the true distance D and the 

brightness distance Dl. Due to the red-shift of the emitted light this becomes [5], 

D L = (l + z)D^, (2.10) 
X 

where 

= sin x , k = 1 

fix) =x ,fc = o 

= sinh x ,k = —1 

Since one observes stars with z S> 0, D and Dl can differ from one another considerably. 



A second possible way of determining distance is to compare the true diameter A of a cosmic source 
with the angle 8 which it subtends at the Earth. In the RW metric, eq.(2.8) implies, 

A D f(y) 

da^j = mm) = —zx (2,11) 

These two examples of how to determine distance show clearly how the space curvature comes into 
astronomical considerations concerning the law of propagation of light. Of course optical methods can 
only be used to determine the distances of objects whose light reaches us. 

2.6 The Friedmann equations 

Having establish the metric (RW), the solution of Einstein's field equations requires a knowledge of the 
form of the energy- momentum tensor, T^ v . To be consistent with the symmetry of the metric we find 
that the simplest realization of such a tensor, T^, is that of a perfect fluid characterized by a time- 
dependent energy density p(t) and pressure p(t). For a perfect fluid symmetry requirements dictate that 
the energy momentum tensor has the form: 

T^ u = (p + v)UfiU v - pg^ u , (2.12) 

where is the fluid four-velocity and for a comoving frame, = (1,0,0,0) ( i.e. a fluid at rest). 
The equation of motion of a particle in a gravitational field is given by 

fjL + r^uV^o (2.13) 

Straightforward but tedious calculations show that the components of the Ricci tensor are [5] 

^oo = -3^, (2.14) 



,R 2R 2 + 2k 
K R + R 2 

and the curvature scalar is 

,R R 2 + k 



U ij = -is + m—)9H, ( 2 - 15 ) 



^ = -6(- + — ). (2.16) 



The spatial 3-dimensional curvature is given by 



3 ^=f (2-17) 



From eqs.(2.12) and (2.14-16), it follows that 



and 

The Einstein's field equations (EFE) are related by the Bianchi identities and only two are independent. 
From eqs.(2.18) and (2.19) we obtain 

| = "f G( P + 3p) (2.20) 

Today R > 0; if in the past p + 3p was always positive, then R was always negative, and thus at finite 
time in the past R must have been equal to zero. This event, referred to as the big bang, is usually taken 
at time t = 0. 

At R = 0, there is a singularity; extrapolating past the singularity is not possible in the framework of 
classical general relativity. 

The constant defined by H = is not a genuine constant but in general varies with time as i -1 and 
hence H -1 sets a time scale for the expansion : R roughly doubles in a Hubble time. 

Equations (2.18) and (2.19) are known as Friedmann's equations. eq.(2.18) can be written as 

im = zmkG- lsa - 1 (2 ' 21) 

where VL = pj p c and p c = 3H 2 /8irG: p c is known as the critical density of the Universe and f2 as the 
density parameter. eqs.(2.18) and (2.20) give 

A = {2q - l)H 2 (2.22) 

where q = —RR/R 2 , a dimensionless parameter, is known as the deceleration parameter which is a 
measure of slowing down the expansion of the Universe, is 

q = ^(l + 3p/p) (2.23) 

Hence qo = ^r(l + 3po/po) 5 where the subscript denotes the present day quantities. The sign of k in 
eqs.(2.21) and (2.22) gives the sign of Q — 1 and q. At the moment observations only yield the bound 
-1< q < 2. 

The parameter k is normalized to the values 

k = +1 — > > 1 closed 

k = -^n = l flat 

k = — 1 ^$7<1 open 

We will consider each of the above cases separately. 



2.7 Successes and shortcomings of the Standard Model 

The standard big-bang model of the Universe had three major successes: 

(i) it predicts that something like Hubble's law must hold for the Universe. 

(ii) it predicts successfully the formation of light atomic nuclei from protons and neutrons a few minutes 
after the big-bang. This prediction gives the correct abundance ratios for 3 He, D, 4 He and 7 Li . 

(iii) it predicts a relic of cosmic background radiation having a black-body spectrum with a temperature 
of 2.75K, today . 

However, certain problems and puzzles remain in the standard model. 

(i) the Universe displays a remarkable degree of large-scale homogeneity. This is evident in the cosmic 
microwave background radiation (CMBR) which is known to be uniform in temperature to about one 
part in 10 4 . 

(ii) a certain amount of inhomogeneity must have existed in the primordial matter to account for the 
clumping of matter in galaxies and cluster of galaxies, etc., that we observe today. Any small inhomo- 
geneity in the primordial matter rapidly grows into a large one with gravitational self-interaction. Thus 
one has to assume a considerable smoothness in the primordial matter to account for the inhomogeneity 
in the scale of galaxies at the present time. The problem becomes acute if one extrapolates to 10~ 43 s 
after the big-bang when one has to assume an unusual smoothness in the initial state of matter. This is 
known as the smoothness problem. 

(iii) the present discrepancy between the observed density of matter and the required value. If O were 
initially equal to unity ( flat Universe) it will stay equal to unity forever. On the other hand, if Q, were 
initially different from unity, its departure from unity will increase with time. 

The present value of Q ranges between 0.1 and 2. For this to be the case the value of $7 would have had 
to be equal to 1 to one part in 10 15 a second or so after the big-bang, which seems an unlikely situation. 
This is called the flatness problem. 

To deal with these problems Alan Guth (1981), proposed a model of the Universe, known as the infla- 
tionary model, which does not differ from the standard model after a fraction of a second or so, but from 
about 1CU 45 to 1CP 30 seconds it has a period of extraordinary expansion, or inflation, during which a 
typical distance (R) increases by a factor of about 10 50 more than the increase that would obtain in the 
standard model. Though inflationary models solve some of the problems of the standard model, they 
throw up problems of their own, which have not all been dealt with in a satisfactory manner [6]. 
The consideration of the Universe in the first second or so calls for a great deal of information from the 
theory of elementary particles, particularly in the inflationary models. This period is referred to as the 
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very early Universe. 

2.8 Matter dominated (MD) Universe 

1. Flat model 

This is the simplest case (k = 0) and is known as the Einstein de Sitter (ES) model. Equations (2.2) 
and (2.12) give 

P + 3|(p + p) = . (2.24) 
The pressure p and the density p are related by the equation of state: 

p = p(p) (2.25) 

which is taken as p = wp for a perfect fluid. 

For dust w = 0, for radiation w= 1/3 and w= -1 for a vacuum dominated Universe. Hence 

p ~ R-W+v) (2.26) 
For the matter dominated epoch p = 0, i.e. w=0, therefore eq.(2.24) gives 

R = R (t/t ) 2 / 3 (2.27) 

where 

Po Rl = P R 3 . (2.28) 

For this case p = p c , $7 = 1. The present value of p c is 2 x 10~ 29 gcm~ 3 . The present age of the Universe, 
to is given by 

to = ^r (2.29) 

which has been found to be less than the range allowed by observations. 

2. closed model 

This is the case k = +1. eqs.(2.20) and (2.18) give 

R 4vr 
R = 3~ 

and 

R 2 _ 8itG 1 
R 2 '^T^R 2 

The age of the Universe, to, is given from the integral 

y/RdR 



Va-R 
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(2.30) 



where a = „ 290 ^ . One has 



+ gO r -l(l-9o) V 2 g0 - Ti . . 

*o = ; —5-777 cos . (2.31) 

Defining the red-shift z as ( z + 1 = eq.(2.30) can be written as 

^ ✓•(l+z)- 1 

1 = H ° Jo [i-n + n 0X -^ (2 ' 32) 

From the above equation we see that the age of the Universe is a decreasing function of f^o : large £lo 
implies faster deceleration, which in turns corresponds to a more rapidly expanding Universe in the past. 
The age of the Universe provides a very powerful constraint to the value of Qq, and to the present energy 
density of the Universe. Independent measurements suggest that 

to = 10 to 20 Gyr 

In terms of z and Qq one can write the age of the Universe as; 

t = !# - 1 (i+*r s/2 , 



for fto = 1 and 



! ftp goz-go+j) 2(l-tt ) 1/2 

* = jff o 777i — — ^-[— cosh( — — — — h 



2(l-O ) 3 /2 L v f2o(l + «) (fi * + l) 1 / 2 fi (l + 2) J ' 

maa: ~~ H (2q„- 1) 3 /2 dL - 2 ffoVn^ 



for < 1. A closed universe has a maximum radius of R mnT , = „ ,„ 2q ° at t 



3. Open model 

This is the case k = -1. eq.(2.18) becomes 



62_,P , ^ a_ 2 g0 



i? 2 = (^ + l), /3 



" ^ #o(l -2g ) 3 / 2 ' 
The age of the Universe is 

_ go r Vl-^gp (1 - go + VI - 2g ) 

t °- J ff (l-2 (7 o) 3 / 2[ go go J " 

Like in the Einstein-de Sitter model, the Universe in this model continues to expand forever. 

2.9 Radiation dominated (RD) Universe 

The solution of eq.(2.24) gives (p = |p) 

p ~ R~ A (2.33) 
10 



and from eq.(2.18), 

R ~ r 1 ' 2 (2.34) 

so that 

p ~ t~ 2 (2.35) 

In this case the curvature term is negligible in comparison with the second term in (2.18). Therefore, k 
= 0, 1, or -1 does not change the physical results. 

The radiation in this epoch is modeled by that of a black body. The density of the black body 
radiation is related to its temperature by 

p r = aT? (2.36) 

where a = 8.418 x 10~ 36 gcm" 3 K~ 4 . 

According to eq.(2.33), the temperature of radiation is inversely proportional to the scale factor of 
the Universe: (This also follows from the fact that the black body radiation retains its spectrum during 
the expansion of the Universe) 

TrOcRT 1 . (2.37) 

When R is very small, that is, in the early Universe, T r can be very high. Equations (2.28) and (2.33) 
give 

f = (f )o(§). (2.38) 

Pm Pm -ft 

This formula shows that the ratio of the radiation and matter densities is not invariant; rather, it 
decreases as the Universe expands. Although the value (p r /p m )o is very small, being only 10~ 3 ; in the 
early Universe, i.e. that is when 

we had 

^»1. 

Pm 

Namely, the radiation was the dominant component of the Universe and its temperature was 

T r = (T r ) (R /R) = T r0 (l + z) . 

For T r Q = 2.7K and z ~ 10 3 , this corresponds to 

T r > 2.7 x 10 3 ~ 3000K. 

This phase of the Universe is called the radiation dominated (RD) Universe. 

Since the Friedmann models are frequently used to interpret cosmological observations, we will now 
derive some of the observable quantities in these models: 
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2.10 The cosmological tests 



2.10.1 The red-shift 

Consider a galaxy Gi at (r, 0, <p) emitting light waves towards us. Let us denote by to the present epoch 
of observation. At what time should a light wave leave G\ in order to arrive at r = at time t = to ?. 
Since a light signal moves in a null geodesic, ds = 0. Then the RW line element (eq.2.6) gives us 

f R dr , 

Since r decreases as t increases along the null geodesic we should take the minus sign in the above 
equation. 

Suppose that the light left Gi at time t\, hence from the above equation 

rU> g]f rri dr 

J h ~W) = Jo VI - kr 2 ' (2 ' 40) 
Thus if we know R{t) and k we know the answer to our question. 

Assume the wave crest was emitted at t\ and t\ + At\ and received at to and to + Ato, respectively. 
Then 

rt +Ato dt ri dr 
iti+At! W)~Jo VI - kr 2 ' (2 ' X) 

If R is a slowly varying function of time, i.e. it effectively remains unchanged over the small intervals 

Ati and Ato we get, by subtracting eqs.(2.40) and (2.41) 

At Ati 



R(to) R(h) 

i.e. 



(2.42) 



The quantity z defines a red-shift (Ai = c5t\,\o = c5t$). The wave length of the light wave increases 
by a factor z in transmission from G\ to us, provided i?(to) > R{t\). Thus Hubble's observations of the 
red-shift is explained if R(t) is an increasing function of time. This red-shift arises due to the passage 
of light through non-Euclidean space time. It does not arises from the Doppler effect, since in our 
coordinate frame all galaxies have constant (r,9,4> ) coordinates. In non-Euclidean space-time it is not 
possible to attach an unambiguous meaning to the relative velocity of two objects separated by a great 
distance. Equation (2.43) may be compared with the gravitational red-shift which is characterized by 
the fact that if light from object B to A is red-shifted, the light from A to B is blue shifted. In the 
present case, if light from A to B is red-shifted, that from B to A will also be red-shifted provided R(t) is 
an increasing function of time. We will therefore refer to the present red-shift as cosmological red-shift. 
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2.10.2 Luminosity distance (Al) 

This is defined by [2] 

D L =nRo(l + z) . (2.44) 

1. Einstein- de Sitter model 
eqs.(2.39) and (2.27) give 

" = l tl R(t) ~ R J tl 
3to 

Rq 



ri = ^[l-(t/to)5] (2.46) 



Upon using eq.(2.29), this becomes 



The luminosity distance becomes 



n = TrV[i-(i + z)- 1/2 ] 



D L = ^-[(l + z)-(l + z) 1 / 2 } 



2. Closed Model 

1 

3. Open Model 



D L = s-afoz + (?o ~ 1)( v 7 ! + 2^o - 1)] (2.47) 



£>L = tt-2 [?o-z + (90 - ^(v 7 ! + 2zq - 1)] (2.48) 
Note that the equations for Dl for the closed and open models are the same. 

2.10.3 The angular size 

We will study how the apparent angular size varies with red-shift in different Friedmann models. We 
will assume that sources of a fixed size d are observed at different red-shifts. Thus a source at (r, 9, <p) 
with red-shift z will subtend at the observer at r = 0, the angle 

(2.49) 

A0i is defined in terms of z and qo. It is interesting to note that A9\ does not steadily decrease as z 
increases. For qo = 1/2 

_ dHo(l + zf/i 

Ml ~ (l + z)l/2-l (2 ' 50) 

A minimum value for A9\ occurs at z = 1.25: 

A6q m i n = 6.75H d. 

The cases (qo > 1/2) are more involved [2]. 
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2.10.4 Source counts 

The number of astronomical sources with red-shifts between z and z + dz is given by (applicable to all 
Friedmann models) [2] 

dN = Wre(t) dr (2.51) 
V 1 — kr 2 

dN = ^ [ ^ + ( g Q-l)(v^p-l)] 2 ^ (2 . 52) 

where n = -j^. 

2.10.5 Particle horizon 

The proper distance to the horizon in a RW space-time is given by 

If R{t) oc t n , then for n < 1, dff(£) is finite and is equal to ^pz^ j i-e. in spite the fact that all physical 
distances approach zero as R — > 0, the expansion of the Universe precludes all but a tiny fraction of the 
volume of the Universe from being in casual contact. This is a vexing feature of the Standard Model. A 
more general expression for the dn [5], is 

dH{t) = H (l + z)Jo [^(l-f^ + ^W 1 - 3 ^] 1 / 2 (2 ' 54) 

From this expression we see that if w < — 1/3, the integral will diverge and the horizon distance will be 
infinite. The present horizon distance is given by 

r 2 k = 0, go = l/2 

° i 7T^ sinh_i yW * = -!,«)< 1/2 
The existence of a finite dn means that the Universe has a particle horizon. 

2.10.6 Event horizon 

A related question to that posed at the beginning of subsection 2.10.1 is whether a light signal sent out 
at the present time (to) reaches all points of the Universe before its end at time Since light travels a 
maximum coordinate distance 

f tE dt 

dEH= L m 

there exists an event horizon if this is smaller than tt or oo: we shall never learn any thing about events 
which at the present time are situated at a distance greater than the above distance [5] 
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Chapter 3 



Thermal History of the Universe 

3.1 Thermodynamics of the Universe 
3.1.1 Introduction 

The first theoretical basis of the evolutionary view is thermodynamics. Applying the requirement of 
thermal equilibrium to the Universe, we have to say that the general tendency of the cosmic evolution 
is for the temperature to be the same all over. It was pointed out by Helmholtz in 1854 that the whole 
Universe will eventually be in a state of uniform temperature and will be falling in the state of thermal 
death (eternal rest). According to the thermodynamics any temperature difference between two systems 
will approach zero. Using thermodynamics, it was proved that for an expanding Universe, even if the 
initial temperature is the same, a temperature difference may still be generated. Consider a spherical 
region R and assume that the matter in R has already reached equilibrium at the beginning. Roughly 
speaking, there are two kinds of matter in the Universe, one is baryonic and the other is radiation. Let 
p m and p r be the mass densities of matter and radiation respectively. Their equations of state are [1] 

Pr = -^Pr (3.1) 
3 

p m ~ m n + -nT (3.2) 

where n is the number density of the nonrelativistics particles and m their rest mass. Here we have 
assumed that each particle has three degrees of freedom. The expansion of the Universe should be 
adiabatic. No exchange with the exterior for the exterior does not exist in the Universe. The expansion 
of the region R which is typical of the Universe, must also be adiabatic. Even though an exterior region 
for R exists, there is no difference between R and its exterior because the Universe is uniform throughout. 
The adiabatic expansion of a system satisfies 

dE = -pdV , (3.3) 
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where E is the total energy of the system. This formula constitutes the basis for the thermodynamics 
of an expanding Universe. We now proceed to investigate the consequences of these equations. 

3.1.2 Radiation under adiabatic expansion 

The total energy density in the region R is 

E r = Vp r , V ~ R 3 . (3.4) 

Substituting eq.(3.4) in (3.3) gives 

d(R 3 p r ) = -p r dR 3 . (3.5) 

Then using eq.(3.1), this becomes 



d(R 3 p r ) = -\p r dR 3 (3.6) 



or 



3 

i.e. 



R 3 dp r + p r dR 3 = -\p r dR 3 (3.7) 



dp r 4 dR 3 



(3.8) 



p r 3 R 3 ' 

The solution of this equation is 

p r oc iT 4 . (3.9) 

According to the thermodynamics of radiation the relation between p r and the temperature of radiation 
T r is 

p r oc T 4 . (3.10) 



From eqs.(3.9) and (3.10) we immediately get 

T r oc R- 1 . 

This result shows that as the Universe expands the radiation temperature falls in inverse ratio to the 
scale factor R , provided only radiation exists. 

3.1.3 Matter under adiabatic expansion 

Assuming the presence of matter only we apply eq.(3.3) to obtain 

d(R 3 Pm ) = - Pm dR 3 , (3.11) 
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where p m is the pressure of matter, 

p m = nT m . (3.12) 

Substituting eqs.(3.12) and (3.2) in (3.11), we then find 

3 

d(R 3 mn) + d(R 3 -nT m ) = -nT m dR 3 , (3.13) 

The total number of particles N = nV is conserved within V, and so n ~ R~ 3 . Thus, eq.(3.13) now 
becomes 

^d(R 3 nT m ) = -nT m dR 3 . (3.14) 



Using N = nV, this becomes 



with the solution 



3 dT m dR 3 



T m oc R- 2 . (3.16) 

This result shows that, as the Universe expands, the matter temperature T m also decreases but in a 
different manner to radiation. 

3.1.4 Generation of temperature differences 

T r and T m vary with R in different ways. Hence during the process of cosmic expansion it is impossible 
for T r and T m to be always equal. Even if initially T r and T m were equal, later on, with expansion, we 
will have T r > T m . Therefore, the cosmic expansion saves the Universe from the final outcome of thermal 
death. According to thermodynamics, a system with thermal equilibrium has the same temperature for 
all the various components, but if radiation and matter always keep the same temperature, how can 
temperature difference ever appear ? 

In a system with complete thermal equilibrium all components must have the same temperature. How- 
ever, a certain time is required to reach equilibrium, interaction between radiation and matter must be 
carried on for a length of time before the two can achieve the same temperature. If the time required 
to achieve uniform temperature is longer than the time scale of the cosmic expansion, then there will 
never be thermal equilibrium between radiation and matter. In this case, it is reasonable to separately 
solve the thermodynamics equations for the two components. Due to cosmic expansion, matter is not in 
a state of complete thermal equilibrium. Radiation and matter are separated from thermal equilibrium 
because there has not been sufficient time for the two components to achieve mutual equilibrium. 
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3.1.5 Equilibrium thermodynamics 

Today the radiation in the Universe is comprised of the 2.75K microwave photons and the cosmic sea of 
1.96K relic neutrinos. Because the early Universe was to a good approximation in thermal equilibrium, 
there should have been other relativistic particles present, with comparable abundance. The number 
density n, energy density p, and pressure of a dilute, weakly interacting gas of particles with g internal 
degrees of freedom, in thermal equilibrium at temperature T, are given by [5] 

n= j£^ I m A ' (3 - 17) 

9 



(27T) 



J E(q)f(q) d 3 q, (3.18) 



and 

> = W)J& nq) *< (319) 

where E 2 = q 2 + m? and / is the Fermi-Dirac (FD) or Bose-Einstein (BE) distribution for species in 
kinetic equilibrium. 

/(g) = (exp^^±l)-\ (3.20) 

where p is the chemical potential of the species, +1 is chosen for FD species and —1 for BE species. It 
follows that 

_ g f°° (E 2 -m 2 y/ 2 E 2 dE 
P ~2n 2 J m exp(V)±l ' 
g f°° (E 2 - m 2 y/ 2 E dE 

71 — TTIo / ,E-„., , '-. ' {6.ZZ) 



2vr 2 J m exp(^) ± 1 ' 
g [<*> {E 2 - m 2 fl 2 dE 

For T 3> m (relativistic limit) and T»/i 



P= £*L exp^±l ' ^ 



f^T 4 , BE 
ImaT 4 FD 



(3.24) 



and 



»= ™„\ (3-25) 



^gT 3 BE 
\£(3)gT 3 FD 

£(3) = 1.20206 is the Riemann zeta function. In the nonrelativistics limit 

mT x3/2„„„ (m-p) 



n = ,(^exp-^ F ^, (3.26) 
p = mn, (3.27) 
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p = nT<^p , (3.28) 
for BE and FD species. For relativistic particles, eqs.(3.21) and (3.23) give 



7T 2 _ 4 



Pr = ^9*T 4 (3.29) 

Pr = \pr = ^2*T 4 (3.30) 
where g* , the total number of effectively massless degrees of freedom (mj «T), is given by 

9*= E ^(^) 4 + ^ E ^) 4 (3.31) 

i=boson i=fermion 

and Tj is the temperature of species i that may not have the same temperature as photons. The factor 
| accounts for the difference in FD and BE statistics. 

3.1.6 Entropy 

Through most of the history of the Universe (in particular the early Universe) the reaction rates of 
particles in the thermal bath, Tj nt , were much greater than the expansion rate, H, and local thermal 
equilibrium should have been maintained. In this case the entropy per comoving volume element remains 
constant. The entropy in a comoving volume provides a very useful fiducial quantity during the expansion 
of the Universe. The second law of thermodynamics implies 

TdS = d(pV) + pdV = d[(p + p)V] - Vdp (3.32) 

The integrability condition [5] 



d 2 S d 2 S 



dTdV dVdT 

relates the energy density and the pressure 



(3.33) 



or 



Equation (3.32) becomes 



T ^=P + P ( 3 - 34 ) 



dp = P 4^dT . (3.35) 



1 ( p + p )v 

dS = -d[(p + p)V] - { -^^dT (3.36) 

dS = d[ ( p+ J^ V + const.] . (3.37) 
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The entropy per comoving volume is (up to an additive constant) 

S = R \ +P) (3.38) 
The general relativistics field equations require 

dS = (3.39) 

This results implies that in thermal equilibrium the entropy, per comoving volume, is constant. It is 
useful to define the entropy density s 

At high temperature the entropy density is dominated by the contribution of the relativistic particles, 
so that to a very good approximationR 

2vr 2 

s = -rr9*sT 3 (3.41) 
45 

9*s= E 9i^f + \ E «(|) 3 - ( 3 - 42 ) 

i=boson i=fermion 

During thermal equilibrium all particle species had a common temperature. Whenever a particle species 
becomes non relativistic and decouples its entropy is transferred to the other relativistic particle species 
still present in the thermal plasma, causing T to decrease slightly less slowly. Massless particles that 
decoupled from the heat bath will not share in the entropy transfer as the temperature drops below the 
mass threshold of a species. Instead, the temperature of the decoupled massless particles scales as R^ 1 . 

3.2 The early radiation era 

The Friedmann equation (2.18) for the radiation era is 

R 2 8vrG . A , 

During this epoch p oc R~ 4 , so that 

i?~t 1/2 . (3.44) 
The energy density of the dominant black-body radiation in this era obeys the equation 

P = \aT\ a=^, (3.45) 



1 s is also proportional to the density of relativistic particles and in particular 
s = 1.8<7* s n 7 , n 7 = photon density 
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where g=2 for photons. Therefore, it follows that 



and that 



The solution of this latter equation is 



R _ T 
R~~T 



2 8irGg*Tr 2 4 



(3.46) 



y T J 3 30 y ' 



giving 



and 



T = ( - )V4 t -l/2 = / ^ \l/4 t -l/2 (- 3 48 N 



T(JT) = ^t-^isec) (3.49) 



To determine which types of particles will be in thermal equilibrium at a given temperature one proceeds 
as follows: we note that if the particle's number density is n, its velocity is v and its reaction cross section 
is a, then the average reaction rate is given by [6] 

r =< nav > (3.51) 

averaged over the thermal distribution. Some relevant reaction are: e + + e~ — > \x + /2, e + + e~ — * v + 9, 
q + q^q + q, £ + l^q + q where q, I stand for quarks and leptons respectively. A given species of 
particles will remain in thermal equilibrium so long as this reaction rate is sufficiently high compared to 
the age of the Universe, i.e. if 

r » H = 1 (3.52) 

Hence initially when n and v are very large we expect all particles to be in equilibrium, but once the 
temperature drops, and the energy of a particular particle falls below the relevant reaction threshold, a 
will vanish and the particle will cease to be in equilibrium. Therefore, when T > H (V < H) the particle 
species couple (decouple) from the thermal bath. 

For reactions mediated by massless gauge bosonsP|[5] V ~ na\v\ ~ a 2 T, during the RD epoch H ~ 

pi 

so that jj ~ - ™ P1 , {\jAira = g =the gauge coupling strength). 

Therefore, for T < a 2 m p i ~ 10 16 GeV, the reactions are occurring rapidly, while for T > a 2 m p i ~ 
10 16 GeV, they are effectively "frozen out." 



2 e.g. e + e+ — > 7 + 7 , e + e+ — > fx + /i H 
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For reactions mediated by massive gauge bosons]]: T ~ na\v\ ~ G\T b and ~ G x m p iT 3 , where Gx ~ 
-%r and mx is the mass of the gauge bosons. Thus for mx > T > G~ 2 / 3 m }^ ~ (m x /100GeV) 4 / 3 MeV, 
the reactions are occurring rapidly, while for T < (rnx/100GeV) 4 / 3 MeV, they are effectively frozen out. 

3.3 Neutrino decoupling 

For T ^> lMeV, neutrons would be in equilibrium with the rest of the Universe via the weak interaction 
processes (e + + e~ — >u + h , ,u + e~^i' + e~ ,e~+p— > v + n , etc.). The weak interaction cross 
section a is given by 

a ~ G 2 F T 2 (3.53) 

where Gf is the Fermi's constant. The number density of neutrinos is oc T 3 . The reaction rate for 
neutrinos (T w i) therefore falls with decreasing temperature as T 5 . Hence eq.(3.51) gives 

% oc T 3 (3.54) 
H 

When eq.(3.52) is not fulfilled, neutrinos decouple (freeze out) from all interactions and begin a free 
expansion. The decoupling of and v T occurs at 3.5MeV whereas the v e decouples at 2.3MeV. At 
decoupling, the neutrinos are still relativistic and thus their energy distribution is given by the Fermi 
distribution and their average temperature equals that of photons. 

As the Universe cools and the energy approaches 0.8MeV, the reactions converting the protons into 
neutrons stop but the neutrons also decay into protons by the beta decay, n — * p + e + u, and therefore 
the ratio of protons to neutrons increases. The mean life of the neutron is 889s and in comparison with 
the age of the Universe which at this time is a few tens of seconds the neutrons are essentially stable. 
The electromagnetic cross sections are oc T~ 2 , and the reaction rate (r em ) is then oc T so that 

^ oc I (3.55) 

Equation (3.52) is satisfied for all temperatures, so, in contrast to the weak interactions, the electromag- 
netic interactions never freeze out. The reaction 

e~ + e + — > 7 + 7 (3.56) 

creates new photons with energy 0.51MeV. This is higher than the ambient photon temperature at that 
time, so the photons population gets reheated. To see this consider eq.(3.41) 

2vr 2 

s = ^9*sT 3 , (3.57) 



3 e.g. e +p^n + v, e +e+^/i ,£ + £^q + q 
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Since entropy is conserved throughout, this is only possible if g* s T 3 remains constant. By applying this 
argument to the situation when positrons and most of the electrons disappear by annihilation below 0.2 
MeV. We denote quantities above this energy by + , and below it by - . Above this energy the particles 
in thermal equilibrium are 7, e~ , e + . Then the entropy is 

s = — aTl . (3.58) 

Below this energy only photons contribute the factor g*s = 2 . Hence, 

T_ = (^) 1 / 3 T+ = 1.40T+ (3.59) 

The number density of neutrinos is related to that of photons as 



n v = \-J[ n ~< ( 3 - 60 ) 



When electrons become slow enough, they are captured into atomic orbits by protons, forming stable 
hydrogen atoms {B.E = 13.6eV). Actually the formation of hydrogen atom occurs at 0.3eV, because the 
released binding energy reheats the remaining electrons, and also because the large amount of entropy 
in the Universe favors free protons and electrons. When this recombination is completed, the photons 
find no more free electrons to scatter against, thus the photons decouple and the Universe becomes 
transparent to radiation [3]. 

3.4 Nucleosynthesis 

We left the story of the decoupling nucleons at the time when the weak interaction ceased and the 
conversion of protons to neutrons stopped because the energy in the thermal bath dropped below 0.8MeV. 
The neutrons and protons were then non-relativistic, so their number densities were each given by the 
Maxwell-Boltzmann distributions. For nuclear species A(Z) with mass number A this is given by 

m A T 3/2 Ai A -m A 
UA = 9A ^^T' eXP ^ f ' ^ ' 

where /j, is the chemical potential of the species. At equilibrium this is related to the chemical potential 
of the proton (jx p ) and neutron (fj, n ) by, 

HA = Z/ip + (A - Z)ji n (3.62) 
By eliminating the exp^^/T) term, one can write [5] 

n A = ^ 3/2 2- i (^) 3(A " 1)/2 <nf Z eMBA/T) . (3.63) 
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Ba = Zrrip + (A — Z)m n — tua , (3.64) 

So the neutron to proton ratio is 

n n /n p = exp(-Q/T) , 

where Q = m n — m p . At energies of the order of (m n — m p ) = 1.293MeV or less, this ratio is dominated 
by the exponential. Thus at T = 0.8MeV one finds that the ratio has dropped from 1 to 1/5. Already 
at a few MeV, nuclear fusion reactions start to build up light elements. For a typical bound nucleus 
A with atomic mass A and number density ua the mass fraction of the various nuclear species A{Z) is 
given by 

Xa = — ^— (3.65) 

n N 

where is the total number density of nucleons. From eq.(3.63) it follows that 

xri^X^X^' 2 'exp(B A /T) (3.66) 

where rj = = 2.68 x 1CT 8 (fish 2 ) , Qb being the baryonic fraction of the critical density. The 
equilibrium abundance of deuterons is given by : 

X d = l6.3(—f 2 v X p X n exp(^) . (3.67) 

Since B^, the deuteron binding energy (= 2.22MeV), is low, X^ is not high enough to start fusion 
reactions leading to 2 H, 3 H and 4 He unless T drops to less than 10 9 K. Although the deuterons are 
formed in very small quantities, they are of crucial importance to the final composition of matter. 
Photons of energy 2.22MeV or more photodisintegrate the deuterons into free protons and neutrons. 
Consider the above process at equilibrium, and define the relative abundance by 

X d 3 n B (2n) 3 V m d } T y ' 

where tun = rriB is the nucleon mass. 

The above equation tells us that as the Universe cools the equilibrium shifts in favor of d over p 
and n at T ~ 10 9 K. We denote the temperature T/v, when the above ratio equals to unity, as the 
nucleosynthesis temperature. 

All evidence suggest that the number density of baryons (nucleons) is today very small. In particular, 
we are able to calculate it up to a multiplicative factor 0_b/i 2 , [3] 

nB = PB_ = ObPc „ 113 x 1Q _5 n B h 2 cm- 3 . (3.69) 
tub m B 
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The parameter fig cannot be very much larger than unity, because that would close the Universe too 
fast. 

Detailed calculations show that deuterons production becomes thermodynamically favorable only 
at energy of magnitude 0.07MeV. Other nuclear reactions also commence at a few MeV. The 3 He is 
produced in the reactions 

d + d^ 3 He + n, 
p + d^ 3 Fe + 7, 

with the binding energy 7.72MeV. This reaction is hampered by the large entropy per nucleon, so it 
becomes thermodynamically favorable only at O.llMeV. The 3 H is produced in the fusion reactions 

n + d^ 3 H + j, 

d + d -> 3 H + p, 
n+ 3 He -^ 3 H +p, 

with the binding energy 8.48MeV. A very stable nucleus is the A He with a very large binding energy 
of 28.3MeV. Once its production is favored by the entropy law, at about 0.28MeV, there are no more 
7-rays sufficiently energetic to photodisintegrate it. A He is mostly produced in the reaction 

d + d^ A He + ~/ (3.70) 

However 3 He and 3 H production is preferred over deuteron fusion, so A He is only produced in a second 
step when these nuclei become abundant. The reactions are then 

n + 3 He -> A He + 7, 

d+ 3 He -> A He + p, 
p+ 3 H^ A He + 7, 
d + 3 H —+ A He + n. 

At energy of magnitude O.lMeV when the temperature is 1.2 x 10 9 K and the time elapsed since the 
big bang is a little over 2 minutes, the beta decay of neutrons already noticeably converts neutrons to 
protons. At this point the n n /n p ratio has reached its final value. 

n n 1 



n p 7 



(3.71) 
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These remaining neutrons have no time to decay before they fuse into deuterons and subsequently into 
4 He. There they stay until today because bound neutrons don't decay. The same number of protons as 
neutrons go into 4 He, and the remaining free protons are the nuclei of future hydrogen atoms. Thus the 
end result of nucleosynthesis taking place between 100 and 700 seconds after the big bang is a Universe 
composed almost entirely of hydrogen and helium nuclei. 

Why not heavier nuclei ? It is an odd circumstance of nature that although there exist stable nuclei 
composed of 1, 2, 3 and 4 nucleons, no nucleus of A = 5 exits, and no stable one with A = 8. In 
between these gaps, there exist the unstable nuclei 6 Li 7 Be and the stable 7 Li. Because of these gaps 
and because 4 He is so strongly bound, nucleosynthesis essentially stops after 4 He production. Only 
small minute quantities of stable nuclei 2 He, 3 He and 7 Li can be produced. The relic abundance of 
the light elements bears quite an important testimony of the Big Bang. The number of 4 He is clearly 
half the number of protons as neutrons go into 4 He. Thus the excess number of protons going into the 
formation of hydrogen is n p — n n . Usually one quotes the ratio of mass in He to total mass in l H and 
4 He, [3] 

4 He 

^itfTTu^ - 25 - (3 - 72) 

At present, the best estimate of I4 from observational data is in the range 0.22 - 0.24 with no directional 
variation. This is a strong support of the big bang hypothesis. The helium mass abundance Y4 depends 
sensitively on several parameters. 

(i) The quantity r] in eq.(3.69) oc Qb- If the ns increases £Ib and i] also increase, and the entropy 
per baryon decreases. Remembering that the large entropy per baryons was the main obstacle to early 
deuteron and helium production, the consequence is that helium production can start earlier. But then 
the neutrons would have had less time to beta decay, so the n n /n p ratio would be larger than ^. It 
follows that more helium will be produced: I4 increases. 

(ii) An increase in the neutron mean life implies a decrease in the reaction rate T w i . Hence a longer 
mean life implies a higher decoupling temperature and an earlier decoupling time. As we have already 
seen, an earlier start of helium production leads to an increase in Y4. 

(iii) The expansion rate H oc ^fgl which in turn depends on the number of neutrino families N u . If there 
were more than 3 neutrino families H would increase . Similarly if the number of neutrinos were very 
different from the number of antineutrinos, contrary to the assumption in standard model, H would 
also increase. From the nucleosynthesis value of r\ one can obtain an estimate of the baryon density 
parameter Ub- Combining eqs.(3.69) and (3.70) one gets 

n B = = c 3.65 x lOV- 2 , (3.73) 

Pc Pc 
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The observational data [3] constrains rj to 

2.8 < 10 10 7] < 4.2. (3.74) 

This limit on rj gives: 

0.010 < Q B h 2 < 0.015. 
Note that h ranges from 0.5 to 0.85 [3] which implies 

0.01 <n B < 0.05. 

Thus we arrive at the very important conclusion that there is too little baryonic matter to close the 
universe. Either the Universe is then indeed open, or there must exist other non-baryonic matter. 

3.5 The decoupling of photons 

The relative abundance of free electrons of number density n e , free protons of number density n p and 
free neutral if-atoms of number density n#(= n B — n e ) in thermal equilibrium at a given temperature 
is determined by 

!^ = ( !|Z:)3/ 2exp( _^ ) (3.75) 

flH 2,71 1 

where B B = 13.59 eV is the hydrogen binding energy. The ionization fraction is given by 

n e 

x = — , 
us 



so that 



x 2 1 ,m e T.o /2 . B H , 



= — ex P(-^)- 



l-x n B 2vr ' " v T 
For f^o^ 2 = 0.1 [2],x = 0.003 , at T = 3000K. Thus by this time most of the free electrons have been 
removed from the cosmological brew, and as a result the Universe becomes effectively transparent to 
radiation. This era signifies the beginning of a new phase when matter and radiation become decoupled. 
This phase continues to the present epoch. During this phase each photon frequency is red-shifted as 
the number density is oc R~ 3 and the temperature is oc i? -1 . A background radiation of temperature 
T ~ 3K therefore means that the red-shift at this epoch when matter decoupled from radiation was 
z ~ 10 3 . The opaqueness of the Universe prevent us from "seeing" directly beyond the red-shift of 
z ~ 10 3 . Thus any evidence of the big bang must come indirectly. The abundance of light nuclei and 
the MBR provide us with the only means of checking the very early history of the Universe. 
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Chapter 4 

Nonstandard Cosmological Models 



4.1 Background 



From eqs.(2.18) and (2.20) we see that if we want a static solution of Einstein's equations, that is, one 
in which R = 0, we must have p + 3p = 0, which is somewhat unphysical, because, if p > then p < 0, 
and if p = then p = 0. Therefore Einstein modified his equations by adding the so-called 'cosmological 
term' to the field equations, as follows 

Knv ~ ^fcg^v - ^9iiu = SnGT^ , (4.1) 

where A is the cosmological constant. Hence eqs.(2.14-16) in (4.1) give 

R 2 k 8vrG A , . 



and 



R 3 vr 1; 3 

The energy conservation law (eq.(2.2)) takes the form 



i? 4vrG . . A ^ , 

(p + 3p) + -. (4.3) 



so that the entropy S, defined in 



dE + pdy = dQ = TdS, (4.5) 



is no longer constant. However, in the standard model the entropy is constant and that was considered 
to be one of the problems of the standard model. This is because the entropy of the Universe, today, is 
unexplainably very huge. Therefore, a variable cosmological term gives rise to an increasing entropy in 
conformity with the second law of thermodynamics. This is true if the cosmological constant decreases 
with expansion. 
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A more general relation between k and f2 is given by 

2A 

fi (l + 3w) = 2q + ^ , (4.6) 

and 

§0(1 + «;) - go - 1 = -J^ . ( 4 - 7 ) 

where w is denned by the equation of the state p = wp, and the subscript denotes the present day 
quantities. Consider a Universe for which k > 0. It then follows for w = that 

«. > 1 - 3^. (4.8) 

Thus we see that the cosmological constant changes the simple relation between k and Q, see eq.(2.21). 
For instance, a flat universe is no longer characterized by = 1. Note that although A is exceedingly 
small, the term may be of the order of unity. The value of f2 = 1 is preferred on theoretical grounds, 
but the observational values are mostly much smaller than 1. 
If we demand that R(t) = i? =cons tant, and p = 0, we get 

" = Sg • k = W a - (4 ' 9) 

Thus A must be positive and therefore k = +1. This is the Einstein's static model. Later on, Einstein 
regretted the addition of this term when he knew about the expansion of the Universe. 

Many dynamical solutions with the cosmological term were studied by Lemaitre and they were 
accordingly known as Lemaitre models. In recent years other motivations have been found for introducing 
a cosmological term. Introducing this term is like introducing a fictious 'fluid' with energy momentum 
tensor given by 

T lv = {p +p')u^u u - p'g^ = ^q9^ , ( 4 - 10 ) 

so that the energy density and pressure of this fluid are given by p' = g^y, p' = — g^j- Hence eq.(4.1) 
can be written as 

- X -n 9ilv = 8irG(T^ + T;j. (4.11) 

4.2 Limits on the cosmological constant 

From eq.(4.4) one gets [32], 

go = ^ - A/3# 2 (4.12) 
\q -n /2\ = \A/3Hi\ 
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Though O an d qo are uncertain one can safely say that — 5 < qo < 5 and < ft < 4, hence |A| = 21Hq. 
By setting Hq = 100kms _1 Mpc~ 1 , this leads to A ~ 10 -54 cm~ 2 . Now in the Newtonian gravitational 
theory if r is the distance of a point mass with respect to the centre of spherically symmetric distribution 
of matter, then the force on this unit mass is given by 

4-7T A 

F = --Gp + -r (4.13) 

where A is the Newtonian form of the cosmological constant: 

(i) A > implies a repulsive force and 

(ii) A < implies an attractive force. 

The matter distribution ceases to be a bound system if F is an outward force. This implies [32] that 

A < AttG P (4.14) 

which gives A ~ 10 _48 cm -2 . From the above we see that A acts at dimensions of order of galaxies and 
has no effect on the solar system. 

4.3 De Sitter model 

Consider a flat RW universe with p(t) = po= c onstant. Equation (4.2) becomes 

R 

— = H = const. (4.15) 
R 

Its solution is 

R(t)ocexpHt (4.16) 

This is known as an inflationary solution. De Sitter obtained this solution for p = p = 0. However, such 
a solution is also possible even without A but with p = const. The de Sitter universe is characterized by 
motion without matter in contrast to the Einstein's universe which is matter without motion. 

One gets a similar form for R(t) as in the steady-state theory of Bondi and Gold (1948) and Hoyle. 
However, unlike in the de Sitter model, which is empty, in the steady state theory there is continuous 
creation of matter due to the so-called c-field [33]. Although we started with k = space-time; inflation 
is not restricted to k = 0. For a Universe with k = +1 one finds the solution [39] 

Ro^H^coshHt , (4.17) 

and for k = — 1 

R(t) <x H' 1 sinh Ht , (4.18) 
which for large t approaches the inflationary solution. 
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4.4 Importance of the cosmological constant 



There is no convincing evidence available for a nonzero value to the cosmological constant. The present 
interest in the flat cosmological constant models has also appeared motivated by two reasons: 

(i) a A term helps to reconcile inflation with observations. 

(ii) with a A term it is possible to obtain, for flat universes, a theoretical age in the observed range, 
even for a high value of Hubble's constant [34]. 

4.4.1 A varying cosmological term cosmology 

The present estimates of the value of A are very small. One way to resolve this dilemma with obser- 
vations is to assume that A is not a 'pure' constant, but rather decreases continuously with cosmic 
expansion. Hence, we can say A is extremely small because the Universe is old. Cosmologists postulate 
a phenomenological law for the decay of this term [26,12]. 

Recently (Waga & Torres), considering the statistics of gravitational lensing, which is a powerful 
tool in constraining models of the Universe, especially those with a A-term, have shown that cosmologies 
with a varying cosmological term give a lower lensing rate. This is due to the fact that in a varying A 
cosmology the distance to an object with red-shift z is smaller than the distance to the same object in 
a constant A model with the same ^o- So, the probability that light coming from the object is affected 
by a foreground galaxy is reduced in a decaying A cosmology [37]. 

Chen and Wu [36] advocated the possibility of particular independence behaviour, i.e. A oc Rr 2 . 
They argued in favor of this behaviour of A from some very general arguments in line with quantum 
cosmology. From dimensional consideration and in the sprit of quantum cosmology, one can always write 
A as Mpj times a dimensionless product of quantities. For an ansatz for the evolutional behaviour of A, 
as in the common practice in quantum cosmology, it is more convenient to use the scale factor R instead 
of the age of the Universe. Supposing that no other parameters are relevant here, the natural ansatz is 
that A varies according to a power law in R. Theoretically, it can be obtained from some simple and 
general assumptions in line with quantum cosmology. Observationally, it is not in conflict with present 
data and may alleviate some problems in reconciling data with the inflationary scenario. 

Therefore, one can write 

where M p i and r p i are the Planck mass and length respectively. With n = 2 in the above equation, both 
Ti and G disappear, since general relativity is a classical theory, and therefore we have A = where 7 
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is a dimensionless number of order unity. The case n = 3 or n = 1 would either lead to too big or too 
small value for the present cosmological constant A(j. 

On the same dimensional grounds Carvalho et al. [42] assumed 

where l p ,t p are respectively the Planck length and time, n is an integer and tn — H~ l is the Hubble 
time. Recalling the general relativity is a classical theory, in order to get rid of ft dependence of n one 
needs to put n = 2. Therefore, A ~ H 2 . 

On the other hand Ozer and Taha [26] proposed a cosmological model with A = jp (i.e. 7 = 1) based 
on the critical energy density assumption. 

Also a cosmological term varying as R~ n , 9/5 < n < 2, was first introduced by Gasperini [35]. 
Consider the Friedmann equation 

rr2 8vrG k A 

For k = 0, A = 3/3iJ 2 , ft = const., the above equation can be written as 

(3 = — — — = const. , 

Pv + P 

where p v = is the vacuum density. This equation is nothing but the Freese et al. ansatz. Hence 
the postulate A = 3(3H 2 is equivalent to the Freese et al. one. One can also write the above Friedmann 
equation as 



where f2 is the density parameter. One can assume that the role of inflation can be described by specifying 

A 
3H 1 



that it drives the Universe to a state of geometrical flatness, corresponding to fi + ^4? = 1. It is useful to 



regard the quantity on the left hand side of the latter equation as £l e ff = 1> with the term regarded 
as the vacuum energy contribution to the density parameter. Using this definition, we see that inflation 
always drives the Universe to f2 e // = 1. 

In this work we will investigate some of these laws and discuss their implications on the evolution of the 
Universe. 

4.5 Variable G models 

4.5.1 Theories in which G varies with time 

Theories of this type were first proposed by Milne, Dirac, and Jordan. Later, Brans and Dicke (BD), 
Hoyle and Narlikar, and Dirac put forward more elaborate theories of this type. A variation of G with 
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time has a considerable effect on the evolution of the Earth and Sun, and on the orbits of the Moon 
and planets. If gravity has changed appreciably over the life time of the Earth, the radius of the Earth 
might have been affected. It has been suggested that the continents all fitted together at one time on 
a much smaller Earth. As the gravitational constant reduced, the Earth expanded to its present size 
and the continents were forced apart. Also a star, like the Sun, in its hydrogen burning phase [3] has a 
luminosity (L) 

L oc G 7 - 7 , (4.19) 

as so would have been brighter in the past if G decreased with time t. The effect of this on life on Earth 
would be enhanced by the fact that the Earth must be moving away from the Sun if G is declining. A 
too fast decline in G would lead to the Sun having already become a red giant. A varying G leads to 
variation in the Moon's distance and period. The orbits of the planets are also modified, and this could 
show up in the radar time-delay experiments [24]. It should be emphasized that some theories in which 
G varies also predict other changes (to preserve energy conservation, e.g.) which can mask the above 
effects. A stronger limit on G follows if 4 .He and 2 H are believed to be synthesized in the big-bang hot 
phase. One then has 

|§| < 10- 12 yr. 1 (4.20) 

The BD cosmology represents the simplest extension of general relativity (GR). In addition to the tensor 
gravitational fields represented by the metric tensor, there is a scalar field (related to the gravitational 
'constant') which is a function of time. BD take A = and seek to satisfy Mach's principle, that local 
inertia properties should be determined by the gravitational field of the rest matter in the Universe, by 
taking 

G- 1 - E 7 (4-21) 

Universe 

The models for k = are particularly simple, since, there 

R oc t q , G oc f 

where 

2(l + «;) 2 

q = , r = 

H 4 + 3u> 4 + 3w 

w being the 'coupling constant' between the scalar field and geometry: w — > oo gives the Einstein- 
de Sitter model. Note that for general w, Gpt 2 = const.. Dirac's 1937 model is obtained by setting 
w = —2/3. However, BD theory makes slightly different predictions from GR for the deflection of light 
by the Sun and for the advance of planets [25]. And if helium is synthesized in the big bang, there is 
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an even stronger limit, w > 100 [24]. Also if G has had greater values in the past than it has now, 
one would expect a small mass density to have the same effects as a larger mass density has later. An 
increasing G would lead to a contraction of the earth. 

4.6 A new formalism of variable G models 

Ever since Dirac first considered the possibility of variable gravitational "constant", G [4], there have 
been numerous modification of general relativity to allow for variable G [5]. These theories have not 
gained wide acceptance. However, recently [6-10] a modification has been proposed treating G and the 
cosmological 'constant' A as non constant coupling scalars. Einstein's equations 

- -Tlg^ = 8itGT^ + Ag^ u , (4.22) 

are considered, where G and A are coupling scalars and the other symbols have their usual meanings 
. The principle of equivalence then demands that only g^ v and not G and A, must enter the equation 
of motions of particles and photons. This approach is appealing because it leaves Einstein's equation 
formally un changed since variation in A is accompanied by a variation in G [6-10], in such a way the 
usual energy conservation law T.^f = 0, holds. If we take the divergence of (5.4) and use the Bianchi 
identities, we get 

G. v T Vil + k, v g vli = (4.23) 

This approach, however, is non covariant and the field equations can not be derived from a Hamiltonian. 
We note that the propagation equations for the scalars are not contained in eqs.(4.24). Despite this 
drawback, there are several advantages of the present approach. This approach could be a limit of a more 
viable fully covariant theory such as five dimensional Kaluza-Klein theory [11]. Various possibilities for 
variable G can be investigated. The problems of the standard model could be solved as in the inflationary 
scenario [12]. The variation of G with time is reasonable, since G couples geometry to matter , and in 
an expanding Universe we expect G = G(t). It is reasonable to assume that the Universe had always the 
Einstein-de Sitter critical density. For unless the Universe had this, it would have diverged, very rapidly 
from it. With the present Hubble constant H p = 5 x 10 -11 yr -1 , the present critical energy density 
P = t&U — 2 x 10~ 47 GeV 4 . The current energy density of the Universe, on the other hand, is between 
10- 47 GeV 4 and 10~ 48 GeV 4 [26]. 

With these assumptions we will show that there are models in which G increases with time and other 
models in which it decreases. 

Recently, the flat Friedmann Robertson Walker (FRW) models have been studied with the present 
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formalism. A number of solutions were presented including de Sitter-type ones relevant to inflation 
[30,31]. 

4.6.1 Bulk viscous solutions 

The role of viscosity in cosmology has been studied by various authors [18,19,29]. It was initially hoped 
that neutrino viscosity could smooth all initial anisotropies and leads to the isotropic Universe that we 
observe today [13-15,22]. The bulk viscosity associated with the Grand Unified Theory (GUT) [6,16,17] 
phase transition can lead to an inflationary universe. It is also known that the introduction of bulk 
viscosity can avoid the big bang singularities [18,19,27]. When viscosity is introduced in a fluid the 
fluid becomes imperfect, i.e. the one in which the pressure, density and velocity vary appreciably over 
distances of the order of a mean free path, or over times of the order of mean free time, or both. For 
such fluids, the kinetic energy is dissipated as heat. For relativistic fluids, the dissipative effects play 
an important role in the history of the early Universe. The energy - momentum tensor of the imperfect 
fluid (T' a p) takes the form [23] 

T' ap = T a(3 + AT a/3 (4.24) 

where (AT Q( g) is regarded as a correction term to the energy momentum-tensor of the perfect fluid T a p. 
In a comoving frame AToo = 0, and in a general frame it satisfies 

u a u^T a() = (4.25) 

The most general form containing viscosity allowed by this condition and the second law of thermody- 
namics is (see Appendix A) 

AT aP = (H^H^Wyt + r]H a ^uJ 7 (4.26) 

where W a p = u a -^ + u/3- a — \g a ^ul n is the shear tensor and H a p = g a p — u a up is the projection tensor on 
a hyperplane normal to u a . r\ and £ are the coefficient of bulk and shear viscosity, respectively. We will 
here consider the bulk viscosity only since the shear viscosity plays no role in a RW model. Therefore 

Hence 

T Lp = (P + P - rjujj)u a up -{p- rjujj)gap 
If we now replace p — rjuJ^ by , the above equation becomes 

T 'af3 = {P + P* )u a Up - p*g af3 

1 The total pressure p* accounts for the isotropic pressure p plus viscous terms can be represented as a polynomial 
in 9 (9 = u7 ) : p* = p — Ylk=x a k9 k where ak are in general functions of p [21,28,41]. 
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Apart from the above replacement, this equation looks the same as that of a perfect fluid. Therefore, 



the introduction of the bulk viscosity does not alter the isotropy and homogeneity of the Universe. 



The field equations with bulk viscosity can be obtained from the general relativistic field equations 



by replacing the pressure term, p, by, p* , where 



p* = p — f]u:, 



(4.27) 



In a Robertson- Walker model, we have uly = 3H, so that p* = p — 3r]H. 
In eq.(4.27) rj is usually taken to have a power law form [19] 



V = VoP 



ri 



(4.28) 



Where t/q > 0, p is the energy density, n is a constant and H is the Hubble constant. In addition to 
the linear dependence of n upon H, as in the above equation, some workers [20,21,28] have considered a 
quadratic dependence upon H, i.e. 



p* = p- 9(H 2 (4.29) 



The C = const models were analyzed by Romero [20] . The bulk viscous models considered so far are 
endowed with particle creation . In the Chapter 7 and 8 we will consider a model with variable G and A 
and bulk viscosity . This combination of G, A and rj has not been considered before. The model turns 
out to have many interesting features. Various models could be reproduced from this model by taking 
particular values of n, where rj = rjop n . 
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Chapter 5 



A Vacuum Decaying Cosmological 
Model 

5.1 Introduction 

The cosmological constant problem [1], that phase transitions in the early Universe would have left the 
cosmological constant larger than the observed upper bound, A > 10- 84 (GeV) 2 [l, 2], by about 120 order 
of magnitude, is a major puzzle of cosmology and particle physics. Considerable efforts were made in 
seeking its solution [1]. One approach attempts to avoid this impasse by allowing A to vary smoothly 
with time so that models in which it was appreciable in the past could be constructed [1]. Examples are 
the models of Ozer and Taha [1], Chen and Wu [4], as well as some of their later generalizations [5,6], 
all of which require A oc R~ 2 , where R is the Robertson- Walker (RW) scale factor. Recently Carvalho, 
Lima and Waga [7] proposed the ansatz 

A = 3/3tf 2 + |j (5.1) 

where (3 and 7 are dimensionless numbers( natural units being used) and H = ^ is Hubble's constant ( 
an overdot denotes time differentiation) . They suggested the (3 -term in Equation (5.1) on the basis of 
simple dimensional arguments consistent with quantum gravity. Special cases of eq.(5.1) when 0=0 are 
: Chen and Wu [4] singular model (7 arbitrary ), Ozer - Taha (OT) [3] nonsingular cosmology (7 = 1), 
and the singularity - free models of Ref. [6] (1/2 < 7 < 1). Carvalho et al [7] have investigated the 
effect of /3-term in eq.(5.1) on the singular model of Chen and Wu [4]. In the present work we study the 
implications of this term for nonsingular cosmologies of Ozer - Taha [3] type. 
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5.2 Field equations 



In a RW universe with a perfect fluid energy - momentum tensor, Einstein's equations with a variable 
A give (a = 3/8vrG) : 

" " = l R ] + W 3 ' (5 ' 2> 

|--*(p + rt = [|]» + ^-|. (5-3) 
In these equations p and p are the cosmic energy density and pressure and k the curvature index. 

5.3 Nonsingular Model 

The introduction of eq.(5.1) and the radiation equation of state V=\p in eqs.(5.2) and (5.3) lead to 

r2 = f 2 ^! + a ° r ~ 2+4 " ' (5 - 4) 

P = ^^R- 2 + «(1 " flAoRT 4 **, (5.5) 

Pv = ^^R- 2 + ^A R^, (5.6) 

Where Aq is a constant and p,, = is the vacuum energy density. Equations (5.4)-(5.6) were analyzed 
by Carvalho et al [7] for the case A > and (3 < 1 corresponding to a singular universe. As pointed by 
them [7] a "natural extension (of their work ) would be to explore different scenarios obtained by making 
Ao < 0," for which the cosmology is nonsingular. One possible line for such an investigation would be to 
consider a nonsingular cosmology of the OT type based on eq.(5.1). This is the theme of what follows. 

5.4 Radiation Universe 

Requiring A < in eq.(5.4) implies 

It is then possible for R to have had an initial minimum non vanishing value at t=0, say. The necessary 
condition for the existence of this minimum in an expanding Universe is R = at t=0. We explore this 
possibility. In order to reach R = at t=0 the exponent of R in eq.(5.4) must be negative, implying 
that (3 < 1/2. Hence from eq.(5.7), 

k < 2 7 . (5.8) 
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Let R=i? at t=0, when R = 0. Clearly R / or else i? 2 (0) = -oo. Thus 

(1 - 2/3) L i?2- 4/3 J < i_2/3' 1 j 

_ g( 7 -/3fc) r (2 7 -fc)(l-/3)^ 
P (l-2/3) J R 2[ ( 7 - f3k)R 2 -^ J ' 1 J 

_ a(7-/3fc) ri /3(2 7 -fc)^ 

P " " (1 - 2/3)# 2 L 1 (7 _ ^#2-4/3 J" 

From eq.(5.10), 

which is the same as eq.(5.9) of Rcf.[6]. Thus the physical condition po > requires. 

k > 7. (5.13) 

Equations (5.8) and (5.12) lead to 

l/2<7<jfe = l, (5.14) 

implying a closed universe. The result (5.14) coincides with eq.(5.13) of Ref. [6]. Classical general 
relativity without or with constant A does not explain the origin of cosmic entropy because classical 
Einstein's equations are purely adiabatic and reversible. In the present cosmology, however, the change 
in the entropy S is related to the temperature T and the change in p v , by [6] 

TdS = -R 3 dp v , (5.15) 

so that from eq.(5.11), 

dS _2a(rr-P) 2/?(2 7 -l)(l-/3)^- 4/3 1 

1 dR~ 1-2(3 [i ( 7 -/?)fl2-40 J ' (5 5j 

in the radiation Universe. Therefore ^ > provided the second term inside the bracket of eq.(5.16) is 

less than unity. Since 2/3 < 1, < (27 — 1 < 1) but > 1) this would be guaranteed, for R > Rq, if we 

simply choose 7 = 1. Henceforth we restrict ourselves to the model with 7 = 1, /3 / 0. We will consider 

two cases separately 

5.5 Nonsingular model: case 1. 

Putting 7 = k = 1 in eqs.(5.9)-(5.12) we obtain 

1 /? 2 " 4/3 

R2 = T^2p [1 -R^ ]<{1 - 2Pr ^ (5 ' 17) 
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p = a(l -P)H 2 , (5.18) 

, - ri ^°" 4 ' 1 (519) 

P " " (1 - 2/?)i?2 [i (1-/3)^-4/3 J' ^- iy i 

po = 0. (5.20) 

these equations extend the 7 = 1, /? = OT model, albeit in a different direction than that of Ref.[6]. 
As in the OT model the initial Universe here also is empty (and cold) but the density parameter 

compared to = 1 in the OT model. The density p attains a maximum 

a a fr. no\ 

Pmax = " — < 77752 ' (5.22) 

2Rl\2(l - 13)}~ 2K o 

at 

R = R max = [2(l-p)}^Ro> Ro- (5.23) 

Note that in Ref.[6] po < aj2R\ and p ma x > cx/ARq. An estimate of Rq and hence an upper bound on 
Pmax can be deduced as follows. From eq.(5.17), 

< R < R 1 , (5.24) 

implying the existence of a natural cosmic acceleration limit in the radiation Universe. Such a maximal 
acceleration, of the order of Planck mass M p i = G -1 / 2 , (h/2n = c = 1), has been discussed before 
[8]. Thus taking R Q 1 ~ M pl yield R ~ G 1 / 2 « 8 x lO-^GeV) -1 = 1.6 x l(T 33 cm, so that p max < 
3 x 10 95 kgm -3 . Finally, for R 3> Ro, R ~ t(linear expansion) and p = p v ~ i~ 2 . The variation p v oc t -2 
is regarded by Berman [9] as more fundamental that expression of p v in terms of R, e.g. p v oc R~ 2 [3-6] 



5.6 Radiation and matter 



Equation (5.2) and (5.3) may be combined to give 



+ !W J»_ (5.25) 



where E = pi? 3 and A may be written as (7 = k = 1): 
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These equations are assumed to hold throughout cosmic evolution. Following Ozer and Taha [3] and 
Ref. [6] we assume the Universe to have evolved through: a very early era, Rq < R < R±, say, of pure 
radiation (discussed in subsection ); a subsequent period of i?i < R < R2 of matter generation; and, 
lastly, for R > R2, an era of radiation and conserved non relativistic matter reaching to the present. 
Except, perhaps, for the matter generation epoch the vacuum is assumed to decay into radiation only. 
Hence with the radiation and matter densities denoted by p r and p m respectively (p r + p m = p), the 
energy, in volume R 3 , of non relativistic matter is E = p m R? = pmpR p = E mp , where the subscript "p" 
designates present-day quantities. We also assume that non relativistic matter has zero pressure so that 
P = Pr = \pr- Chen and Wu [4] and also Carvalho et al [7] assume that the vacuum decays into non 
relativistic matter in the matter-dominated (MD) universe. But Freese et al [10] have demonstrated 
that such a scenario is not, apparently, favored by observations. Under our assumptions eqs.(5.25) and 
(5.26) yield, when R>R 2 , 

dp r 4(1 - 0) _ 3(3 Pmp Rl 2a(l - 0) 

~dR + ~^~ Pr ~ R± + ~ R3 • ^- 27 ) 



This has the solution 



where 



^M + 7 ^ [1 + ^, , ( , 28) 



(1-40)R 3 (1-20)R 21 R?~W 



(l-20)R 2 p p m Prp 30 

1 } a(l-0) [ Pmp l-4/3 J 
3d 

= (1 - 20) [ Prp - -^—p^l - 0) Pvp - 0( Prp + Pmp )]-\ (5.29) 

Hence from eq.(5.26), 

_ /Wg q(l - 0) 0uR^ 

Pv (1-40)R 3 + (l-20)Ri [ + (1-8)R 2 ~^ U 



and by eq.(5.2), 



-1 



From eq.(5.31), 

n = l-0<l, (5.32) 
as before (see eq.(5.21)). Lastly, the present value of the deceleration parameter q = —RR/R 2 is 
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For (3 = eqs.(5.28)-(5.33) reduce to the corresponding OT results. Denote by R eq the value of R at 
t = t eq , the time when radiation and matter were balanced in equilibrium. Let E r = p r R 3 , where p r is 
given by eq.(5.28), be the radiation energy. Then 

E r {Req) = E mp (5.34) 

(5.35) 



E r (R) > E mp when R < R eq 
E r (R) < E mp when R > R cq 

The condition that E r was decreasing as R approached R eq would then imply that 



^2-4/3 

(l-4/3)u,>-^. (5.36) 
Rp 

Hence either 

to > and (3 < 1/4, or, lo < and 1/4 < (3 < 1/2. (5.37) 
We discuss each case separately . 



1. u > 0, f3 < 1/4 



From eq.(5.29) we have 



Hence 



3/3 < Prp_ 

1-4(3 p mp ~ ' 



<^ = §rp . (5.38) 

Pmp 

(3 < 5 rp . (5.39) 



There are three subcases: 

(a) (3 = 0: This corresponds to the OT [3] model. 

(b) (3 > : Assuming 5 rp <C 1, i.e. the present Universe is MD, we have < (3 <C 1. Also from eq.(5.30) 
dp v /dR < 0, so that entropy is continuously generated. From eq.(5.28), 

6r - p m " 1 - 4/3 + (1 - 2(3) Pmp Rl X [1 + H*-V J - (5 - 40) 
In the early radiation and matter universe S r 3> 1. Using this condition in eq.(5.40) and noting from 
(5.38) that 

3/3(1 - 4/?)" 1 < <5 rp « 1 one has 

(l-2/W?3 

+ i? 2 -^ ^ a(l-/3)i? • 1 j 



But by condition (5.36) 
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for R eq /R > 1 and < /3 < S rp < I. Hence from eqs.(5.31), (5.41) and (5.42) , 

R ~ (1 - 2P)R*-*> ~ (1 - J) ' (5 ' 43j 
when R eq /R 3> 1. Consider now p„ of eq.(5.30). Assuming, plausibly, that for R <C R eq the vacuum 
energy was much more important than the energy of non relativistic matter, i.e. p v /p m = S v S> 1, we 
have 

A^- 4 ^ (i - 

(1 - ^i? 2 -^ ^ a(l - " 1 j 



But by (5.36), 



(1 - 0)R 2 -^ (1 - (3)(l - 4(3)R 2 ~ 4 P ' 
For R <C i?eq, and provided (3 is not negligibly small ( much smaller than 5 rp ) , the RHS of (45) will be 
considerably larger that unity. ( If (3 <C S rp <C 1 , the /3-term in A may be dropped, which is subcase 
(a)). Then eqs.(5.30), (5.45), (5.44) and (5.43) imply 

^ a/3u,B*-*P ^ (3p r 
Pv ~ (1 - 2/3)^-4/3 - (TT^) • t 5 - 4 ^ 

The relation between p v and p r in this result is identical to eq.(5.7) in the paper of Freese et al. [10] 
where a parameter x = p v /(Pr + Pv) replaces (3. But whereas letting x — ► in Ref. [10] produces 
standard cosmology, one is not entitled to the limit (3 — ► in eq.(5.46) here because the derivation of 
this equation assumed that /3 is not vanishingly small. 

(Taking naively (3 — > in eq.(5.46) one concludes erroneously that p v = in the OT model). Integration 
of eq.(5.43), taking approximately t = ti ~ when i? = i?2 ~ 0, gives 

i? = [ 2 ( 1 -^ 1 ^" 2/? ] (l/2(l-/3) t l/(2(l-ffl (547) 

Hence from eq.(5.46), 

pv " 4(1 - m 2 o^W) pr ' ( } 

which is the same as eq.(5.8) of Freese et al. [10]. As noted by them in eq.(5.48) has been suggested 

by various authors [11], with (3 presumably dependent on the particular model of vacuum decay. In the 
present model < (3 < S p . 

(c) (3 = — | (3 |< : Then one has from condition (5.36), 

\R\, , 7^2+41/31 

1 _ \P\ ujR P f5 49) 

(1+ | p \)R 2 +\P\ l °- 4yj 

< l _ I P I -Keg 



(l+\(3\)(l + A\P\)R 2 +m ' 
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Thus provided | (3 | is not vanishingly small both sides of this inequality will be negative when R ~ 
R2,Req/R S> 1. (If | (3 | is vanishingly small subcase (a) is retrieved). It then follows from eq.(5.30) 
that p v < and ^jf > for R ~ R 2 , implying a decreasing entropy in the early radiation and matter 
universe. Hence we exclude (3 < 

Observe that eq.(5.31) implies R < if oj > and (3 < 1/4. On the other hand from eq.(5.24) R > in 
the radiation universe. This reversal of the sign of R is suggestive of an intermediate phase transition 
period separating the pure radiation and the radiation and matter eras. 

2. u = - I w |< , 1/4 < (3 < 1/2. 

Consider eq.(5.31). Combining the second and third terms on its RHS and using condition (5.36) give 

1 r l^l^" 4/3 1 , , 

h ' (5.50) 



(1 - 2/3) 1 R 2 ~^ 
1 R 2 -^ 



(1-2/3 1 (4/3- l)i? 2 - 4 /3 J ' 
For R2 < R < i?eq t ne expression to the right of the inequality, and so also that to its left, is negative. 
This leads to R 2 < and p < in eqs.(5.31) and (5.31), which is physically inadmissible. Thus we 
conclude from this subsection that to > and < [3 < 5 rp , where S rp is present ratio of radiation-to- 
matter energy density. 

5.7 Phase transition 

From eqs.(5.25) and (5.26) one has 

dE 

— = 2a(l - (3)+ 3(3R 2 p + -3(1 - (3)R 2 p , (5.51) 
leading, on using E r = p r R 3 and eq.(28), to 

"R-2 



rR'2 

3 / R 2 [{1 - (3)p - (3p]dR = 2a(l - (3){R 2 - R ) - E 2 

J Rn 



- wi MR a-^P 1 a(l-/3)(l-4/3) J R 2 FH 

- -2a(l - /3)i?o - (1 _ ^ + , (5.52) 



where 



(1 - A(3)R 2 2 



But by condition (5.36), 



^2-4/3 

< 1 ^ — ™ < 0. (5.54) 

(1 - 4/3) 2 ^- 4 / 3 
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Hence the integral in eq.(5.52) in negative implying that 



^ < 3/3(1 - (3)- 1 < 1, (/3 < 1/4) , 

or by condition (5.38), 

^ < 3/3(1 - I3Y 1 < 3/3(1 - 4/3)- 1 < 5 rp , (5.55) 

during at least part of the period (i?i,i?2). The condition (5.55) is clearly incompatible with the pure 
radiation equation of state. The condition (5.55) is incompatible with the equation ^ = 1 — ^g- of 
radiation and matter. For if one were to suppose that it is consistent with this equation one would have 
1 — ^ff- < 5 rp , which does not admit simultaneously the assumptions of a present MD universe (5 rp <C 1) 
and an early radiation - dominated universe ( E m /E 

Thus we interpret the era (i?i, R2) as a phase transition period associated with the creation of rest - 
mass and the appearance of decelerated cosmic expansion. During this period the pressure may become 
negative. (If (3 = negative pressure must occur.) 



5.8 Baryon-to-photon ratio 



If the radiation produced by vacuum has a Planckian thermal distribution then its temperature T will 
be related to its density p r by 



7T 2 _ 4 



Pr = ^5e//T 4 . (5.56) 

where g e ff is the effective number of spin degrees of freedom. Then in the early (i? <C R e q) radiation 
and matter universe we have from eq.(5.56) and (5.48), 

T(t) * [ 2^r^y^ ]1/4t " /2 ■ (5 - 57) 

which together with eq.(5.47) imply 

T ~ R 13 - 1 . (5.58) 

Equations (5.57) and (5.58) are identical to the corresponding results of Freese et al. [10] . If the photon 
Planckian spectrum is to be maintained the energy per photon G 7 must red-shift like the temperature 
so that G 7 ~ R 13 ^ 1 . Then eqs.(5.48) and (5.47) imply that the photon number density 
riy oc pry/ e y ~ R 3 ^' 1 ^. On the other hand the baryon number density ng ~ i?~ 3 . Hence the baryon-to- 
photon ratio 

^=M^ J R-3^^ T 3/J/(1-/J) > R^R eq . (5.59) 



This is eq.(12) of Freese et al. . 
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5.9 Cosmic nucleosynthesis 



Primordial synthesis of the light elements nuclei crucially depends on the temperature-time relation and 
the expansion rate of the early Universe. In the present model the temperature-time relation is given by 
eq.(5.57). The expansion rate, from eq.(5.43), is 

coinciding with the rate deduced in Ref.[10]. Primordial nucleosynthesis depends also on rj (eq.(5.59) 
here). Since eqs.(5.57)-(5.60) coincide with the corresponding relations of Freese et al. [10] they would 
lead to these authors nucleosynthesis constraints on j3, namely /3 < 0.1 . But the present models require 
< (3 < S rp . Hence it is consistent with the observed element abundance from primordial nucleosynthesis 
provided 5 rp < 0.1 - a condition that readily obtained in a present MD universe. 

5.10 Entropy generation 

The decreases of r] with temperature displayed by eq.(5.59) is a consequence of the production of the 
entropy by the vacuum. Requiring that r] at nucleosynthesis falls in the range 10~ 10 < rj < 10~ 9 , and 
that at grand unification i]gut < 10~ 4 , Freese et al. [10] obtain (3 < 0.12, which is comparable to the 
constraint in the preceding section. 

5.11 Deceleration parameter 

Define 

t = *r P ~ YZJp (5-61) 
By condition (5.38) £ > . Hence from eq.(5.33) , 

0<^<[^ (1 -^- 4 ^ ][l + %^f]--<l-^. 
More precisely since < (3 < 6 rp <S 1, then £ ^ 1 and < q p < 1/2. 

5.12 Age of the Universe 

By Equation (5.29) and (5.61) 



(3mp£, 



(1-2/3)^ (i -/?)(! + u) ' 



(5.63) 
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and from eq.(5.31), 

(1-/3) 

Let R = uR p . For all u > U2 = R2/R P eq. (5.31) can be written as 



ff 2 = Hmp^-rvrpl _ 



,2 H P r e , (1-/?) , ^ 



U " (l-<5 rp ) [ l + ^ + (1-4/3)*, + (1 + lj)v?-W ] ' (5 ' 65) 
Taking approximately t = t2 ~ when -u = U2 ~ we may write for the age of the Universe 

It is readily verifiable that t p (w; /3), /3 < 1/4, is a decreasing function of u) so that 

H p t p {u- 0) < H p t p (0; (3) = (1 + 6) 1 / 2 J, (5.67) 

where 

j= /V/ 2 [i + M~ 1/2 ^ 

JO 

= (1 + b) 1 / 2 yT+6-y/b 

and 6 = ^rgrp ^ 1- Thus J — ► (1 + 6) -1 / 2 or H p t p (0; (3) 1 so that one has the upper bound 
tp(u;P) < Hp 1 » 10 42 (GeV)- 1 « 2 x 10 10 yr 

5.13 Classical cosmological tests 

Returning to eq.(5.31) using eqs.(5.63) to replace p mp by A and introducing the red-shift z(l+z = R p /R) 
we have 

R=(l- 2/3)- 1 / 2 [a(l + z)[l + A^(l - /^(l - 4/3) (1 + z) 1 " 4 ' 3 ] + l] 1 / 2 , (5.68) 
where by eq.(5.63), (5.64) and (5.61), 



a = w(l-/?)(l-4/3)- 1 A = 



;i-/3)(l-2/3)fl 2 ff 2 
(l-4/3)(l + 5 rp ) 



(1 - 2(3)RlHl / x 

(1 + V P ^ 5 - 69 ) 



Now A 1 = < £ <C 1, so that for optical and radio cosmic source with z < 5 ( z < 1 if normal galaxies 
and not quasars are considered) [15] ), A _1 (l - /3) _1 (1 - ^(3){l + z) l ~ 4 ^ <C 1. Also a = (1 + w)^ 1 > 1 
, i.e. RpHp 3> 1. This means that although /3 <C 1, the /3-term in A p is still important, vis-a-vis the 
7-term. Thus eq.(5.68) may be approximated simply by 

R*i v^(l-2/3)- 1/2 (l + z) 1/2 . (5.70) 
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The RW metric can be written as 

ds 2 = dt 2 - R 2 [d X 2 + r 2 { X )[d9 2 + sin 2 9d(f> 2 }} (5.71) 

where r(\) = sinx when k = 1. For light propagation along a radial ( dd = d<p = 0) null geodesic from 
a distant cosmic source at (x, 0) to us at x = eqs.(5.70) and (5.71) give 

(5.72) 

which upon integration yields (x(z = 0) = 0): 

X= 2(1 "^ )1/2 [l-(l + ^)- 1/2 ]- (5-73) 
V d 



Because a S> l,sinx ~ X- Hence 



2(1 + 6) 1/2 



e = ^ P r(x) « ^ ' [1 - (1 + z)- 1 / 2 ] 

= (l + b) 1/2 HES, (5.74) 

where £es denotes the product R p r(x) in the k = Einstein-de Sitter model. This result is important 
for the following cosmological tests: 

5.13.1 Magnitude versus red-shift relation 

The apparent and absolute bolometric magnitudes m and M of a galaxy of red-shift z are related by [16] 

m- M = fj, = 5log 10 D + 25, (5.75) 

where (i is the distance modulus and D = £(1 + z) is the luminosity distance of the galaxy, measured in 
Mpc. From eq.(5.74), 

M~A^s + -log(l + 6). (5.76) 

5.13.2 Galactic diameters 

A galaxy of linear dimension has angular diameter 

d6 = (1 + (5-77) 

Hence by eq.(5.74), 

dfl ss <W ES (1 + 6)~ 1 / 2 . (5.78) 
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5.13.3 Number counts of sources 



The number of uniformly distributed optical and radio sources in the volume element dV , 

dV = R 3 r 2 ( x )dttdx, dQ = sin 6 d9d(/), (5.79) 

is [17]: 

dN = nR 3 r 2 (x)dndx = n p Rlr 2 {x)dVLdx (5.80) 
where n(t) oc p m (t) is the number density and nR? oc p m R 3 = p mp R p so that nR 3 = n p R 3 . Noting that 

n p _ p rnp _ (1-/3) _ (1 - 4/3) 



n* s p m p ES (l + S rp ) (1 + 6)' 
and using eqs.(5.80), (5.74), (5.72) and (5.69) we obtain 

(1 - 4/3)n^ 5 £| 5 (l + VfPdzdSl 



(5.81) 



dN « ' p ," OK -— . (5.82) 

H p (l + z)3/2 V 7 

On the other hand in the MD universe of the Einstein-de Sitter model 

d N ES = n^&sRfdSldx = n p S ^sdzdn 

H p {l + z)2 

Hence 

dN » (1 - 4/3)(l + b) l / 2 dN ES (5.84) 

Equation (5.76), (5.78) and (5.84) are the model's predictions for the classical cosmological test. Because, 
b, [3 <C 1 these results are approximately the same as in the ES model. 

5.14 Nonsingular model: case 2. 

Here we will take the case for which 7 = 1-/3 and k = 1, so that po / 0. Equations (5.9), (5.10) and 
(5.11) become 

R 2 = l- A 2 ~^ , (5.85) 
R 

P=^[l-(l-/?)(§^- 4/3 ], (5-86) 
Pv = ^[1 - /3(^)] , (5.87) 

And for physical reason we will take (5 > 0. The rate of change of entropy at temperature T is given by 
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= 2a[l- 2/3(1 -/?)(§) 2 - 4 ^] (5.89) 



The equation shows that ^ > for all R > Rq and thereby solving the entropy problem of the standard 
model . 

The density parameter Q, defined by 
Q = pj p c , p c = aH 2 becomes 

= 1-/3+4 ( 5 - 9 °) 
R z 

The maximum of p is attained at 

R = R mx = [2(l-f3) 1 ^-W)} (5.91) 

Pmx = p{R mx ) = ^^ ) ^- (5.92) 

Using eq.(5.9) one can estimate the value of Rq since R = < Ro 1 and that signal a limit for a 

cosmic acceleration in the early Universe. Such as limit is of the order of Planck mass (M p i = G -1 / 2 ). 
This Universe accelerates less rapidly that the one with 7 = 1. Thus starting with some non zero initial 
density allows the Universe to be less accelerating than with zero one; and with a less maximum density 
than the one with 7 = 1. 

5.15 Radiation and matter 

In the wake of pure radiation era, the rest mass is generated during the period, say, R\ < R < R2 [6]. 
For R > i?2, after creation of the rest mass, the matter energy (E m ), E m = p m R 3 stayed constant so 
that E m = PmR? = E mp = PmpRp, where "p" denotes present-day quantities. Thus for the present case 
with A given by eq.(5.1) and 7=1-/3, k = 1 one gets 

Pr = ^|ir 3 + aR~\l + , (5.93) 

where u + l = a^ 1 p mp Rl(5 p - and 8 P = p rp /p mp - 

The radiation energy E r , E r = p r R 3 is given by 

E r = 3 J^ + aR[l + u(^r^], (5.94) 

and 

~ T^p ~ T^I/3~ + T^p { li ) ■ (5 - 96) 
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The acceleration parameter, q, 

, = =^ ( 5 . 97) 

The present value is 

l/2 + [u,a(l-4P)/(l-f3)]p mp Rj 

ip i + [( i _ p + w)(1 _ W( i - /3)] PmpjR 2 - • 

Denote by i? eg the value of R at t = t eg , the time when radiation and matter were balanced in equilibrium 

E r (R e q) = E m p- 

The condition that E r was decreasing as R approached R eq would then imply that 

(1 _ 4/3)0, > (^) 2 ^ . (5.99) 
Thus either to > and /3 < 1/4 as w < and 1/4 < (3 < 1/2. 

However, the later case implies that R 2 < 0, p < for R2 < R < R eq - In Ref.[3] one approximately has 
p r R 4 constant or OT constant for R < R eq in the radiation and matter universe. In the present case for 
R <C R eq the first two terms in eq.(5.95) contribute negligibly to p r and therefore 

* w^ 2 ' 1 " • (5100) 

and 

R 2 = T^(f) 2 ~ 4/3 • ( 5 - 102 ) 

This leads to 

Pv = JZTpPr, R<^Re q . (5.103) 

This is the same relation as postulated by Freese et al. [10] which suggests that in the early phase of 
radiation and matter p r and p v red-shift at the same rate. The solution of eqs.(5.102) and (5.103) is 

R ~ tV2(i-« , (5.104) 

and 

Pv = —pPr « 4{1 _ p)2 f2 ■ (5-105) 
valid in the early radiation and matter (,i.e. R <C R eq ) epochs. 

This behaviour of p v ~ t~ 2 has been noted by many workers. Though the model of 7 = 1, k = 1 and 
7 = 1-/3 has different features at the early Universe, they give the same results in the wake of early 
radiation and matter phases of the Universe that are identical to Freese et al. 
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5.16 Phase transition 



We consider, finally, the matter generation period R\ < R < R 2 . For R > R 2 eq.(5.95) implies R > 
if oj > and (3 < 1/4 in contrast with R > for R < R\ as previously noted. Thus the appearance of 
rest-mass ushers in decelerated expansion of the Universe. 
Equations (5.25) and (5.26) with 7=1-/3, k = 1 

— = 2a(l - 2/3) + 3(3R 2 p - 3(1 - (3)pR 2 (5.106) 
ait 



and using eq.(5.94) one gets 
IRo 

(E 2 = E m {R p ) + E r (R 2 )) 



rR-2 

3/ R 2 [(l-f3)p-f3p}dR = 2a(l-2(3){R 2 -R ) + E -E 2 , (5.107) 



where 



2a(l-2/3)(i? 2 -i? )+/3a«o-p mp ^- r ^P mp ^-ai? 2 [l+^(^) 2 - 4 ^ , (5.108) 
= -a(2 - 5/3)i? - Y-^fe P fip + a(l - 4/3)F(l^) , (5.109) 

as long as (3 < 1/4 the above integral is negative, implying that 3p/p <C 1 for some values of R in the 
interval (Ri,R 2 ). During this period the Universe undergone a phase transition, an era separating the 
pure radiation and the radiation and matter epochs. It was noted by Freese et al. that the constraint 
that an early radiation epoch be followed by a matter dominated era requires that x = (3 < \ in both 
matter and radiation epochs. Nucleosynthesis in this model proceeds the same as in the previous one. 

5.17 Discussion and concluding remarks 

Our aim has been to discuss a nonsingular Robertson- Walker cosmological model with a varying cosmo- 
logical constant A = ^Rr 2 + 3(3H 2 . This form of A is due to Carvalho et al. [7]. In the considered 
model the Universe was initially (t=0) empty and had a minimum scale factor Rq > 0. Subsequently it 
evolved through consecutive phase of pure radiation, matter generation and radiation and matter. The 
absence of the initial singularity requires 2(3 < 1 and \ < 7 < k = 1. In the very early pure radiation 
era R > 0. For R S> Ro in this era R ~ t, compared to R ~ in the standard model. On the other 
hand generation of entropy throughout the pure radiation era requires [see eq.(16)] 2/3 ^ 2 J~^^~^ < 1- 
To satisfy this condition we have chosen, for simplicity, 7 = 1 and have, henceforth, confined our selves 
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to this case . The model is thus an extension of the Ozer-Taha model(/3 = 0,7 = 1) albeit in a different 
direction than that followed by Ref. [6]. General physical consideration and an increasing entropy in 
the radiation and matter period place a stronger restriction on (3, viz. < (3 < 5 rp , where 5 rp is the 
present ratio of radiation to matter energy densities. This constraint is quite stringent if one assumes, 
as we have done here, that the present Universe is matter -dominated so that (S^ < 1. Under these 
conditions the pure radiation and the radiation matter eras are separated by a phase transition period 
that reverses the sign of R from R > to R < 0, ushering in decelerated expansion. But the precise 
phase transition equation of state is not known. In the early radiation and matter universe the model 
is virtually indistinguishable from the flat-space decaying-vacuum singular cosmology of Freese et al. 
[10]. The condition < (3 < 5 rp means that the nucleosynthesis constraint of Freese et al, (3 < 0.1, 
is easily satisfied for 5 rp <C 1. This provides an interesting connection between the present and early 
Universe. The Ozer-Taha model and the model of Freese et al. are very dissimilar, as already noted 
by the latter authors [10]. It is therefore remarkable that an extension of the one converges to the 
other. Freese et al. make the basic postulate that the vacuum and radiation energy densities red-shift 
at the same rate, for large R. This feature emerges here in the very early radiation and matter era, 
as a mere by-product of the approach. In the cosmology of Carvalho et al. the Freese et al. scenario 
corresponds to 7 = k = 0, (3 7^ 0. Here also (3^0, but 7 = k = 1. As to the age of the Universe and 
the deceleration parameter, our model predicts t p < H~ x and < q p < 1/2. On the other hand, in the 
singular cosmology of Carvalho et al. (3 is a free parameter that can be adjusted to produce t p > H~ l 
and q p < 0. We have also examined the implications of the model for the classical cosmological tests. 
The results approach Einstein-de Sitter model predictions. In the present work the density parameter 
$7 = 1 — (3 = constant, before and after the phase transition. Hence 1 — 5 rp < $7 < 1 although this implies 
the absence of flatness fine-tuning problem it would apparently require the existence of dark matter. 
Finally, the proposed cosmology predicts a continuously expanding Universe. 
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Chapter 6 

A Nonsingular Viscous Cosmological 
Model 

6.1 Introduction 

In one approach for solving the cosmological constant problem [1, 2] cosmologies with decaying vacuum 
energy are introduced ( [3,4] and references therein). But in most of the proposed models dissipative 
effects that were probably very important in the early Universe are not considered. On the other hand 
the influence of bulk viscosity on cosmic evolution in standard Friedmann models has been discussed in 
many works( for an excellent review see [5]). In these works the coefficient of bulk viscosity r], say, is 
usually taken to have a power law dependence on the cosmic energy density p, viz. rj oc p n , n constant. 
Recently Calvao, de Oliveira, Pavdn and Salim [6] generalized the model of Chen and Wu [7], in which the 
vacuum energy decreases with cosmic expansion as the inverse square of the Robertson- Walker (RW) 
scale factor, by appending to the energy-momentum tensor a bulk viscosity term with rj oc p. The 
resulting field equations in this case are found to have no analytical solutions in general and methods 
of qualitative analysis of differential equations are used to deduce the cosmological implications of the 
model. 

The present paper also attempts to incorporate bulk viscosity in a decaying vacuum cosmology. But 
unlike Ref. [6] we take n = 1/2 in the dependence of rj on p, i.e. rj oc p 1 / 2 . For this case the solutions of 
the field equations for flat cosmological models with vanishing cosmological constant A appears to have 
a special status: 

They are the only solutions with structural stability [8] . In addition, Beesham [9] has shown that the 
n = 1/2 models are equivalent to the flat variable-A cosmologies of Berman [10]. Because of these 
peculiarities of the models with rj oc p 1 / 2 we limit our self to this case. 

We also assume a Universe with exactly the critical density. Consequently the vacuum energy decays 
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as in the Ozer - Taha model [11]. The critical density assumption, which is indicated by inflation, is 
theoretically appealing because the apparent closeness of the present and critical densities would else be 
difficult to understand on other than anthropic grounds. 

Our model which is nonsingular and closed is introduced in Sec. 5. 2 . Several implications of this 
model turn out to be equivalent to the corresponding consequences in previously proposed non viscous 
cosmologies [12 — 14]. 

6.2 The Model 

In a homogeneous isotropic universe with a perfect fluid energy-momentum tensor, Einstein's equations 
with a variable cosmological 'constant' A give (k = 3/87rG)[ll]: 

r2 -i A k 

d(pR 3 ) 9 k odA , 

where p and p are the cosmic energy density and pressure, respectively, and k the curvature index. To 
obtain the field equations with bulk viscosity we replace p by the effective pressure [15] 

p* = p - 3r]H, (6.3) 

where H = is the Hubble's constant and r\ is the coefficient of bulk viscosity. 
As remarked in the introduction we set p equal to the critical density. 

p = Pc = kH 2 , (6.4) 

throughout cosmic evolution. Consequently eq.(6.10) gives the Ozer- Taha [11] result 

A=| (6.5) 

The form A oc R~ 2 has also been suggested by Chen and Wu [7] on the basis of dimensional arguments 
consistent with quantum gravity. Pavdn [16], on the other hand, has used the Landau-Lifshitz fluctua- 
tion theory to study the physical consistency of several A variations and has concluded that A should 
vary as R~ 2 . The detailed observational consequences of this behaviour have been examined in Refs. 
[7, 11, 13, 17, 18]. As mentioned in the introduction we choose rj oc p 1 / 2 . Specifically we write 

11 = 20j^ pl/2 ' m "° ' C ° nSt " (6 ' 6) 
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Together with eq.(6.13) and the usual equation of state 
p=(7 — l)p,0<7<2 this leads to 

- & " 1 " (6 ' 7> 

where Mpi = G -1 / 2 is the Planck mass. Equation (6.11), with p — > p* , and eqs.(6.16) and (6.14) give 

dp 3 r? n 2k/c 

^ + i? [7 -M^ ]/?= l^' (6 - 8) 

which has the solution 

k~ 1 P R 2 = R 2 = 37 _ 2 2 fc _ + Air^ +2+4 ^ (6.9) 

where /3 = 4 ^° ; and ^4 is a constant. In particular for the radiation- (7 = 4/3) dominated (RD) phase 
of the Universe 

R 2 = -A- + AR- 2+ ^. (6.10) 

1 

Carvalho, Lima and Waga [12] have proposed a model Universe based on the cosmological constant 
ansatz 

3a 

A = 3(3H 2 + , (6.11) 

where a and /? are dimensionless constants of order unity. The detailed cosmological implications of this 
postulate have been explored in [12,3, 14]. In particular a nonsingular scenario based on eq.(6.19) was 
studied in [14]. 

It follows from eq.(6.20) that in the RD universe [12] 

R 2 = 2 ^ + A R~ 2+ ^, (6.12) 

and 

P = n^fjr^ 2 + «(! " /3)A^- 4+4/3 - (6.13) 
For a = k eqs.(6.21) and (6.19) are formally equivalent. It is therefore interesting to ask to what extent 
their equivalence is reflected in the corresponding cosmologies. To answer this question we consider a 
nonsingular scenario as in Ref. [14]. 

A Universe with a non vanishing minimum scale factor Rq at t = arises from eq.(6.19) if A < 0, f3 < 1/2 
and k = 1 (closed universe). Then 

P= lP = R?{1 - 2/3) [ ~ (6 - M) 
with po = and p maximum (= 2 (i-/3)Ri ) a * 

R = R mx = [2(1-/3)] ^i? - 
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An estimate of Rq may be obtained as follows. From eq.(6.19) < R < Rq 1 so that Rq 1 represents 
a cosmic acceleration limit. Such a limit of the order of Mpi has been suggested before [19]. Taking 
R 1 ~ Mpi yields Ro ~ lCT 33 cm. 

We next investigate the implications of the model in the cosmic phase of noninteracting radiation of 
density p r and pressureless non relativistic matter of density p m (p r + p m = p). 

We assume that the vacuum couples to radiation only [11, 13, 14, 18]. Hence p = p r /3 and the matter en- 
ergy E m = p m R 3 = pmpR? p = E mp , where subscript "p" denotes present-day quantities. Then eqs.(6.11), 
(6.14) and (6.16) give (k=l) 

dpr , 4(1-/3) 2k 4(3 

m + ^- pr = l# + m Emp (6 - 15) 



which has the solution 

4(3E mp 



+ 



rr (1-4/3)^3 (1-2(3)R 2 
where 



uRl-W 
1 -i E 



(6.16) 



(1 +u) = (1 - 2[3) K - 1 p mp R 2 p [5 p - j^Lp], (6.17) 

5 P = y 2 - being the present ratio of radiation-to-matter energy densities. It follows from eq.(6.24) and 
conservation of non relativistic matter that 

R ^ K - lpR ^0^ + _^l_ + (1 _ m - K (6 . 18) 

Denote R at t = t eq , the time when radiation and matter were equal, by R eq , i.e. E r (R eq ) = E m (R eq ) = 
E mp , where E r (R) = p r R 3 is the radiation energy. Then the condition that E r was decreasing as R 
approached R eq leads to 

^2-4/3 

(l-4/?)u;>-|^3. (6.19) 
Rp 

Implying that either uj > and /3<l/4orc<;<0 and (3 > 1/4. But the latter case would, by virtue of 
eq.(6.26), imply that R 2 < for R < R eq and is therefore excluded . 

Equation (6.27) is formally identical with eq.(6.24) in Ref. [14]. But whereas (3 there determines the 
coefficient of H 2 term in the cosmological constant ansatz, here it is a measure of bulk viscosity . 
The condition (3 < 1/4 implies, from eq.(6.26), that R < compared to R > in the pure radiation 
era. Assuming that a period of rest-mass generation R\ < R < R2 , say, separated the pure radiation 
and the matter and radiation epochs one can readily show, in exactly the same way as in [14], that this 
period is a phase transition era during which the pressure becomes small or negative. 
The conditions u > 0, (3 < 1/4 and eq.(6.25) lead to 4(3 < 5 p . Hence if the present Universe is matter- 
dominated (MD) < 4S p /3 <C 1 , or, from eq.(6.16), that the present bulk viscosity contribution to 
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the effective pressure p* is at most of the order of the radiation pressure. To examine nucleosynthesis in 
the present cosmology we proceed as follows. Noting that p mp Rp = p r%eq R\ q where p rfiq is the radiation 
density at equilibrium we obtain from eqs.(6.26), (6.24) and (6.27) the approximate relation 

k " * (1 - 2P)R 2 -^ * K ' lprR2 ' (6 ' 20) 
valid for R <C R eq , (3 <C 1 and p m <C p r . 

Equation (6.28) is approximately the standard (k = 0) model's cosmic expansion rate. Counting 
three neutrino species one therefore has for the neutron-to-proton "freeze-out " temperature Tpw{= 
Tp/10 10 K) « 1.11, the same as in standard cosmology. Also from eq.(6.28) p r /p r ~ —4(1 — (3)R/R « 
—4R/R leading to the standard temperature-time relation. 
From eqs.(6.28) and (6.25), using p r = ^g e ffT 4 , one has 

a-i^ = [I^iMi-a ^1/(1-/3) 



■ l + u; 1/4 r 



(6.21) 



where 

Aral-4/W" 1 . (6.22) 
Denoting the baryon number density by ns we then have 

-b««b P [^] 3/4 [|] 3 (6-23) 

for T < T^, the nucleosynthesis temperature, and where ub p = Km^QsHp, with being the nucleon 
mass and £1 b the present baryonic fraction of the critical density. The equation determining T/v (inGeV 
) therefore becomes 

|^ + \ In T N + ln(ft B /i 2 ) - 14.47 = 0, (6.24) 
In 2 

where Bd{= 2.23MeV) is the deuteron binding energy, h is the normalized Hubble constant (0.4 < h < 
0.8 [20]) and 

^B = [^] 3/4 n B >n B , (6.25) 

with the equality sign holding for the non viscous (A — ► 1) standard (w — ► oo) case. 

Apart from the replacement £Ib — ► ^b, eq.(6.32) coincides with the corresponding standard model 

relation [21]. Requiring to satisfy the standard model constraints on i.e. [22] 

0.015 < ^b < 0.070 (6.26) 
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gives TV (standard). With Tp and T/y as in standard cosmology and the nuclear physics aspects of 
nucleosynthesis unaltered the standard nucleosynthesis scenario is reproduced. The condition (6.34) 
implies that 

nB <A 3 / 4 <1, (6.27) 



0.070 

irrespective of the size of the contribution of vacuum energy. In the presence of viscosity (A < 1) the 
constraint (6.35) implies a stronger upper bound on the baryonic density than required by the standard 
model. 

6.3 Concluding remarks 

We have presented a homogeneous isotropic cosmological model with decaying vacuum energy and bulk 
viscosity. Our aim has been to generalize an earlier work by Calvao et al. [6] in which the bulk viscosity 
coefficient depends linearly on the cosmic energy density p. In the present model it varies as p 1 / 2 . 
Our model is nonsingular, closed and, in several aspects, equivalent to previously proposed non viscous 
variable-A models [12 — 14]. This type of equivalence of viscous and non viscous variable-A cosmologies 
has been noted and discussed by Beesham [9]. 

A noteworthy feature of the present work is that primordial nucleosynthesis proceeds as in standard 
cosmology provided the Universe today is matter dominated. This condition implies that the present 
bulk viscosity contribution to the cosmic pressure is, at most, of the order of the radiation pressure. 
Another consequence of nucleosynthesis is that the upper limit on the baryon density is lower than the 
standard value. This strengthens the case for nonbaryonic dark matter. 
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Chapter 7 

A Viscous Universe with Variable G and 



The role of viscosity in cosmology has been studied by several authors [1-4]. The bulk viscosity associated 
with grand unified theory phase transition can lead to inflationary universe scenario. It was well known 
that in an early stage of the Universe when neutrino decoupling occurred, the matter behaves like viscous 
fluid [22]. The coefficient of viscosity is known to decrease as the Universe expands. Beesham [20] 
studied a universe consisting of a cosmological constant (A ~ t~ 2 ) and bulk viscosity. He showed that 
the Berman model could be viscous model for n = 1/2 . 

More recently Abdel Rahman considered a model in which the gravitational constant, G, varies with 
time but energy was conserved [11]. Other models have been considered in the literature by Sistero and 
Kalligas et al . In the present work we will investigate the effect of viscosity in a universe where G and 
A vary in such a way that energy is conserved . 




7.1 



Introduction 



7.2 



The Model 



In a Robertson Walker universe 



dr 2 = dt 2 - 



R\t)[ 



dr 2 



+ r 2 (d6 2 + sm6 2 # 2 



)] 



(7.1) 



1 — kr 2 



where k is the curvature index. 



Einstein's field equations with time dependent cosmological and gravitational "constants' 




(7.2) 
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and the perfect fluid energy momentum tensor 



Ttxv = {p + P)u^u v - pg^y , (7.3) 

yield the two independent equations 

3 | = -4, G (3p + p-4), (7.4) 

3 f = 8 ' G <" + 4>-§- < 7 ' 5 > 

Elimination of R gives 

3(p + p)R = -(^p + p + Jj^)R . (7.6) 
The conservation of energy and momentum yields 

3(p + p) = -R^. (7.7) 

The effect of bulk viscosity in the field equation is to replace p by p — 3r]H, where rj is the viscosity 
coefficient. It follows immediately that 

9 v HR=(^p + ^)R, (7.8) 

and 

p + 3H{p + p) = . (7.9) 

eq.(8) can be written as 

9 4=§" + 8^- <7io) 

where prime denotes derivative w.r.t. scale factor R while dot is the derivative w.r.t to cosmic time t. 
In what follows we will consider a flat universe, k = 
Equation (7.5) and (7.17) lead to 

8irGp = 3(1 - /3)H 2 , (7.11) 

and the equation of state 

p=( n -l)p (7.12) 

in eq.(7.8) and (7.9) lead to 

p = AR-^ , (7.13) 

where A is a constant. 

9r/| = - // + , (7.14) 
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or 

#\ *J_ 1 _ 9 V0 A^ 37 n+3 7 -l 

H 2 2(1 + P)RH 2(1 + /?) 
where we have taken the viscosity coefficient to have the power law 

7] = 7]op n , 7]q > , n const. 

and the ansatz [30] 

A = 3(3H 2 , (3 const.. 

The solution of eq.(7.15) is obtained as follows 
Let y = jj and a = 37 ^~^ , Therefore 

d = -9(1-^)^^-1 n+37 _ a _ 1 

dir 2 

9(1 -^4_^_ 37n+37 _ a 
^ 2(3 7 n-3 7 + a) 



and finally 

eq.(7.11) and (7.21) give 
where 



9(l-^A-- (37n _ 37 + o)J; _3 T „ H .3 7| 



= 2(3 7 n - 3 7 + a) 37TO _ 37 
9(1 - (3) m A"-i 

= 3^(1-/3) ^3^-1) 
8ttA 



D= 2(2n-l-/3) 



9(1 - (3) m A n - 



eq.(7.21) gives 



i i 



R(t) = [3D 7 (1 - n)] 37d-n) i 3 7 (l-n) . 

Hence eq.(7.13) and (7.22) become 

p(t) = A't^z , 
G(t) = B't^r , 

and 

r){t) = A t^ , 

where 

A' = A[3m 7 (l - n)]T^, 
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B' = §^[3m7(l -n)}~ , and A = A' n r] . 
The Hubble parameter is 

where < n < 1. This condition on n rules out some models with n > 1 [3]. The cosmological 
constant becomes 

A = W^?? <7 ' 28) 

This law of variation of A is thought to be fundamental [20]. The vacuum energy density (p v ) is given 

by 

- = • <7 ' 29) 

and from eq.(7.11) and (7.17) we obtain 

Pv = jzrpP- ( 7 - 3 °) 
For an expanding Universe, i.e. H > 0, we must have D > 0. This implies that 

2n - /3 - 1 > 0, 

or 

(3<2n-l. (7.31) 

Whether G increases or decreases depends on the value of n. For n > 1/2 G increases with time and for 
n < 1/2 G decreases with time and for n = 1/2 G remains constant. The condition (7.31) now gives 
A > 0, G increases with time 
and 

A < 0, G decreases with time or remains constant. 

To solve the persisting age problem of the ES model we must have t H > 2/3. This result require that 
G is an increasing function of time. Recently, Massa (1995) proposed a model in which G increases with 
time. An increasing G would cause the Planck length lp = \fG to be an increasing function of time, 
and the quantum fluctuations on the metric would be vanishingly small in the very early Universe. A 
fully classical description of the Universe for all t > would be possible. This is perhaps one of the 
reasons to consider the increasing G model [6]. 

Following Freese et al. one can put a stringent constrain on the value of (3 . The parameter x of 
Freese et al. is equivalent to (3 and since x < 0.1 this implies (3 < 0.1 for the nucleosynthesis constraints 
to hold. The deceleration parameter is given by 

RR 

q = ~lv ' 
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q = 3 7 - 3-fn - 1 . (7.32) 

This shows that the deceleration parameter is constant. The constant deceleration models have consid- 
ered by Berman and Som [9,15]. Equation (7.27) can be written as 



and for the present phase "p " 



It is evident that negative q p would increase of the present age of the Universe. 
From eq.(7.25) we obtain 



G 2n - 1 1 



G 1-n t ' 



(7.35) 



and the present value is 



,G, 2n-ll 2n-l, . „ ,„ . 



A power law dependence of G was obtained by Kalligas et al. [19], it has been shown to lead naturally 
to A ~ t~ 2 . Unlike the model of Abdel Rahman and Beesham, this model shows a constant G does 
not imply constant A. We see that the quantity Gp satisfies the condition for a Machian cosmological 
solution, i.e. Gp ~ H 2 , (see [25]). This also follows from the model of Kalligas et al . 
The relationship between our model and that due to Kalligas et al. is manifested in the following re- 
placement 

1 + n K 

n = 7: 

2 + n K 

and 

2 + n K 

where 

tjk : n due to Kalligas et al . This furnishes the resemblance. Hence Kalligas et al. model is equivalent 
to a viscous model . 

7.3 The horizon problem 

The horizon distance,i.e. the size of the causally connected region, is given by 

n , \ A* dt' 

dH = m LWY 
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We would like to have 37 — 3jn = 1, so 

37-I 



n = 



3 7 



Note that 37 — 3771 > implies n < 1. 

In what follows we will discuss some classes of models. 

7.4 Model with n=l 

Equation (7.7) becomes 



„ _ 9(1-^-^. 3 
a 

y = = const , (7.39) 



a 
or 

Hence 

R(t)=FexpH t , F const., 

which is an inflationary solution. 

Such a solution has been obtained by several authors [2,5,13]. Here the density is not constant but has 
the following variation 

p = AF~^ exp-37# i- 

Such a solution was obtained by Berman and Som for the Brans-Dicke theory for the scalar field <p where 
4> = 1/G [26]. 

In the present case, however, G is not constant during this epoch, viz. 

G(t) = Mexp37# i, 

where 

3H$F*r(l-0) 
8irA 
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7.5 Model with n= 1/2, 7=1 

Equations(7.23)-(7.25) become 

R(t) = Ft§, F const. 
/9 = J 4o^ 2 5 ^0 const. 

and 

G = const. 

The Hubble parameter is H(t) = ft -1 . This is the flat FRW universe result. 
The deceleration parameter is 

q = -3n + 2 = 1/2 

Since several authors claim that the age of the Universe computed from the FRW flat model tends to 
be smaller than the range given by observation, 0.6 < H p t p < 1.4, our model could give a better value 
for any departure from n = 1/2 . 

However, it was found that only n = 1/2 solution are structurally stable [21]. 

It was shown by Beesham that Berman solution (a power law for R) is a viscous solution with n = 1/2. 
The relationship between our model and Berman's [9] is 

m = 37(1 — n) . 

The value of m in our case is not put by hand, but emerges naturally from the dependence of the viscosity 
on the energy density(?7 ~ p n ) in a given epoch. This solution seems more elegant. 

7.6 Model with n=0 , 7=1 

Equations(7.23)-(7.25) give 

R(t) = F't 1 / 3 , 
p = A't-\ 

and 

G(t) = Bt' 1 , 

where F', A' and B are constants. This is a model of constant bulk viscosity. It resembles the Brans-Dickc 
model [16]. We see that 

g = _l = _ 3i y-i, (7.41) 
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(§)p = -ZH p l . (7.42) 

This solution was obtained by Berman [15] for the Bertolami equation for the present phase. 
Note that in GR, k = 

' = sbr <7 - 43) 

Whether our result is acceptable or not depends upon the value we measure for (£) p for the present 
phase. 

This also resembles the Dirac no creation model. For this class of solutions q = 2. 



7.7 Model with n=2/3, 7 = 1 

The scale factor is given by 

R(t) = Rot Rq const., 

and 

P = A't~ 3 , 
G = Bt, 

where A' and B are constants. 

This linear variation of G has been found by Berman [15] for the Bertolami solution for the Brans- 
Dicke theory (BD) with a time varying cosmological constant for the present phase. For this model 
q = 0. 

7.8 Model with n=3/4, 7 = 4/3 

The scale factor is given by 

R(t) = Ft, 

and 

p = AT\ 
G = Bt 2 , 

where F, A 1 and B are constants. 

This solution was obtained by Berman [15] for the Bertolami theory for the radiation era. He also found 
that T oc R -1 , preserving Stefan's law. It was also found by Abdel Rahman that a variable G and A 
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model lead to a similar result for the radiation universe [11]. In his model he considered A ~ R~ 2 . For 
this class of models q = 0. 

More recently, (1995) Massa has considered a model which support an increasing G constant. In 
his work his considered a "maximal power hypothesis (MPH) in the Einstein-Cartan theory of gravity. 
Equivalence of his work and ours requires n > 1/2. 

7.9 Model with n=l/2 , 7 = 4/3 

For this model 

R(t) = Ft 1 ' 2 , 
p = A!t~ 2 , 

and 

G = B = const., 

where F, A! and B are constants. 

This special value for n gives a constant G in both radiation and matter epochs. This is equivalent to a 
FRW flat universe. For this class q = 1. 



7.10 Model with 77 = r] H 



Using eq.(7.13) in eq.(7.14) we obtain 



Let y = 1/H . This becomes 



H' 1 37(1-/?) 1 9(1 ~ P)m pa-y-i (7Ad , 
JP + 2R JP = 2A R (7 ' 44) 



dR R A 



(7.46) 



9(1-/3)^ 37 

^ ^(2a-3 7 ) ' 1 J 



and hence 

H 2 = NR^, N = ^^pl. 

9(1 - /?)% 

Substituting this in eq.(7.11) and using eq.(7.13) yields 

^ 3iV 

G = = const.. 

8ttA 
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(7.48) 



Hence 

R(t) = N^t&. 

This reduces to the flat FRW model with constant G. This is equivalent to the solution with n = 1/2. 
Therefore the assumption 77 ~ H is equivalent to 77 ~ p 1 ! 2 [12]. 

7.11 Model with n = 2/3, 7 = 2 

The scale factor is given by 

R( t ) = Ft 1/2 , 

and 

P (t) = A'r 3 , 

G(t) = Bt , 
7](t) = r^r 2 , 

where F, A' and B are constants. 

This result is obtainable from Berman [2] if we let A = y^U, B 

7.12 Model with n = 1/2, 7 = 2 

The scale factor is given by 

R = Ft 1/3 

and 

P = A't~ 2 , 
G = B = const., 

where F and A' are constants. 
This is the solution for the BD theory for the present phase, as shown by Berman and Som (1990). For 
this model q = 2. This solution also found by Beesham for Bianchi type I models for n = (where 
7 = n + 2, i.e. p ~ i~ 7 ). Barrow showed that p ~ i~ 7 dominates the viscous term for all fluids with 
1 < 7 < 2 [3]. 



= \ and m = 2. For this model q = 1 
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7.13 Model with n= 1/2, 7 = 2/3 



The scale factor is given by 

R = Ft , 

and 

p = A't-" , 
G = const., 

where F and A' are constants. 

These are the solutions obtained by Pimentel [14] for the scalar field of the second-self creation theory 
proposed by Barber, assuming a power law of the scalar field and the expansion factor. The resemblance 
is evident if we put 

_ np + 3(7 - 1) 
67 

np : n due to Pimentel. When n = 1/2, the present case, np = 3. For this class of model q = 0. There 
is no horizon problem associated with this solution. 



7.14 Model with n=l/2 7=1/3 

The scale factor is given by 

R = Ft 2 , 

and 

P = A't~ 2 , 
G = const., 

where F and A' are constants. 

This is the wall- like matter. For this solution q = —1/2. This solution has been obtained by Berman for 
the radiation universe, i.e. a wall-like matter behaves the same as radiation in a viscous universe. 

7.15 Model with n = 0, 7=1/3 

The scale factor is given by 



R = Ft, 
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and 

p = A't-\ 

G = Br 1 , 

where F, A' and B are constants. 

This solution which solves for the power law is also Machian i.e. Gp ~ H 2 , ( see [25]). Such solution 
has been noted by Berman and Som for the constant deceleration type with m/0 [25]. This solution 
corresponds to the case m = 1. In this case we see that the viscosity is constant, i.e. 77 = 770- This 
solution is a wall-like matter solution. This model is free of the horizon problem. This solution has been 
obtained by Pimentel [24] for the solution of the Brans-Dicke theory with a constant bulk viscosity for 
k 7^ solution. He has shown that these solutions satisfy the Machian condition and the second Dirac 
hypothesis. Singh and Devi [23] studied cosmological solutions in Brans-Dicke theory involving particle 
creation and obtained a similar solution for k = 0. Some other solutions are 

7.16 The Pimentel solution for the scalar tetradic theory 
A 

Case I 



This solution of Pimentel [4] is equivalent to our solution provided 

n P = (p P - 2 + 3 7 ) 

and (5p = 2 — 67(1 — n), where the subscript "P" is the Pimentel value. Therefore this solution is a 

(2np-/3p+2) 

viscous solution. The viscosity (7?) varies as t ( 2 -/3p) . The condition (3p / 2 is equivalent ton / 1. 
Note that 

7-2fl- 1 



UJTlp 



or 

2 



to 



3 7 (2n-l)( 7 -2) 
Case II 



This solution is equivalent to our solution provided we make the following substitution 

3 7 (1 - n ) = 1 

or 

2n-l . . 

np = = 37 — 2 

1 — n 
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and therefore 

2(n - 1) 



ijj = 



(2n-l)( 7 -2) 
The viscosity coefficient in this case varies as 

^t-(i+np)( / g p = 0). 

The viscosity term in this case ~ exp — ( "^* )t. 

7.17 Berman solution 

Berman studied a constant deceleration model [2]. eqs. (7.14) and (7.15) of Berman are equivalent to 
eq.(7.24) provided 

(3 = 2n - 1 < 0, 

A 1 

127r 7 2 (l-n)' 

and 

B 2n-l 



3 7 2 (l-n) 2 ' 

We conclude that Berman solution is equivalent to a bulk viscous model with variable G and A. The 
viscosity term here varies as n ~ t~^ 1+ ^\ 

More recently Johri and Desikan (1994) have considered cosmological models in Brans-Dicke theory with 
constant deceleration parameter. Their solution for a flat universe, [their eqs.(65)-(67)] are equivalent to 
our solution, i.e. eqs.(7.23)-(7.25) for the replacement of 
(3 = 37 — 3771 — 1 and a = 3 7 (1 — In). 

7.18 Cosmological expansion in the presence of quadratic 
bulk viscosity (£) 

This term appears as 3£H 2 in the pressure term. 
Let us consider £ = const, case. It follows that 



R G r 8ttG 
Using eq.(7.4) we obtain 

, 37(1-/3) _ 9(1-/3)% 3^3^! 2 _ 3C(1 ~ P) d3t-1 Q-3 = q 

2R 2 2 
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This equation admits a power law solution of the form 



H = aR m , a const. 



Substituting this in the above equation , 

37 
2 



we get m = — t? and n = 1/2 



Hence 



or 



3 7 

H = aR-~. 



„ ,37a 2 JL 
R = t*t 

This is the familiar FRW flat universe solution. If we take a general power law for (, i.e. ( ~ p r for 
some r, it follows that only r = 1 is possible. This case has been studied by Wolf [5]. He showed 
that a constant p leads to the inflationary solution. This model is similar to the one considered before 
{v = Vop)- We see for all these models one has the Machian solution Gp ~ H 2 . 

7.19 Brans-Dicke solution 

In the Brans-Dicke theory [29] the scalar field is related to the gravitational constant G as 4> oc 1/G. 
This theory is equivalent to a bulk viscous solution with 

n = — , (7.49) 

6(l+a;) ' v ' 

in which case the viscosity coefficient varies as 

Zuj+2 

7] ~ t 3 ^+ 4 . (7.50) 

As in the Brans-Dicke theory when u — > 00 the theory approximates to Einstein de Sitter, in the present 
case we have n — > 1/2. 
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Chapter 8 

A flat Viscous Universe with Increasing 

G 

8.1 Introduction 

In most variable G models [7,8] G is a decreasing function of time. But the possibility of increasing G 
has only been considered recently [2]. Massa [5] proposed a model in scale of increasing G depending 
on a "maximal power hypothesis (MPH)." Recently [1], we have considered a cosmological model with 
variable G and A and bulk viscosity. Various solutions are listed and all of them satisfy Mach's condition 
[4]- 

8.2 The Model 

In a Robertson Walker universe 

dr 2 = dt 2 - R 2 (t)[ 2 + r 2 (d6 2 + sin0 2 # 2 )] (8.1) 

1 KT 

where k is the curvature index. 

Einstein's field equations with time dependent cosmological and gravitational "constants" 

~ \g^R = SirGT^ + A 9tlu (8.2) 
and the perfect fluid energy momentum tensor 

Tpv = (p + p)u^u u - pg^ v (8.3) 

yield the two independent equations 

3 | = -4, G (3 P + ,-^), (8.4) 
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4 = 8 * g <" + 4>-§- < 8 - 5 > 

Elimination of R gives 

3(p + p)fl = -(|p + p + (8.6) 
The conservation of energy and momentum yields 

3( P + P ) = -R^. (8.7) 

The effect of bulk viscosity in the field equation is to replace p by p — 3rjH, where rj is the viscosity 
coefficient. It follows immediately that 

9r ) HR=(^p + -^)R (8.8) 

and 

p + 3H(p+p) = 0. (8.9) 

Equation (8.8) can be written as 

"4 = §"+8Sr <8 - 10) 

where prime denotes derivative w.r.t. scale factor R and dot is the derivative w.r.t. cosmic time t. In 
what follow we will consider the flat universe, k = 0. Equation (8.5) becomes 

3H 2 = 8nGp + A. (8.11) 

We take the Chen and Wu ansatz for A [3] 

A = a const., (8.12) 
R 

and the viscosity to have the form 

V = VoP n , % > 0, n const.. (8.13) 

Using the equation of the state p = (7 — l)p in eq.(8.9) we obtain 

p = AR-^, A const.. (8.14) 

Substituting eqs.(8.12), (8.13) and (8.14) in (8.10) yields 

G 1 - %QA n - x R-^ n +^- x HG - ^-i? 37 ~ 3 = 0. (8.15) 

4mA 

This equation admit a power law solution of the form 

G = FR m , F,m const.. (8.16) 
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Inserting this in the above equation we get 

H = C i? 3 ^™- 1 ) (8.17) 

and the condition, 

m = 3 7 - 2, (8.18) 

where C = (m - i^V^o^™ -1 )- From eqs.(8.16) and (8.18) we see that in the radiation epoch 
(7 = 4/3) G (x R 2 and in the matter epoch G oc R. As long as 7 > 2/3, G continues to increase. 

This result agrees with that obtained by Abdel Rahman [2] for the critical density model. Hence his 
model is equivalent to a viscous model. Now eq.(8.17) gives 

R = [3 7 (1 - n )C][l/37(l-n)] t [l/37(l-n)] ( 8 . 19 ) 

Hence eqs.(8.16), (8.14) and (8.12) become 

G = G t [(37 " 2)/37(1 ^ n)1 , (8.20) 

p = P0 i [ - 1/(1 - n)1 , (8.21) 

A = Aotl" 2 / 3 ^ 1 -™)] , (8.22) 

and 

H = w^a ■ (8 - 23) 

where Go,po and Ao are constants. For an expanding Universe, i.e. H > 0, we must have C > 0. This 
implies a < ^irFA(3j — 2), and for a > 0, 7 > 2/3. Since some authors think that the variation of 
A oc t~ 2 is essential, we impose the condition 

3 7 (1 - n) = 1. (8.24) 

Hence eqs.(8.19) and (8.20) become 

R = R t, R const. (8.25) 

G = G t [(2n ~ 1)/(1 - n)1 , G const.. (8.26) 

This solution has been obtained by Ref.l for a different ansatz of A. In the next section we will consider 
two cases separately. 
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8.3 I. Radiation dominated (RD) universe 

This is characterized by p = \p. The condition (8.24) gives n = 3/4. Hence, eqs.(8.21), (8.20) and (8.13) 
become 

/)ocr 4 , (8.27) 

G oc t 2 , (8.28) 

and 

r\ oc t~ 3 . (8.29) 

We see that G increases with time. In Ref.2, for t 3> Rq the gravitational constant G oc t 2 . Though the 
model of Ref.2 is nonsingular and closed yet it evolves towards this singular viscous model. One also 
notices that T oc R^ 1 a result that is expected to hold in this era. An increasing G is recently supported 
by Massa [6]. 

8.4 II. Matter dominated (MD) universe 

This is a dust filled universe (7 = 1). The condition (8.24) gives n = 2/3. Hence, eqs.(8.21), (8.20) and 
(8.13) become 

p oc r 3 , (8.30) 
Goct , (8.31) 

and 

r\ oc t~ 2 . (8.32) 

We see that G increases linearly with time. This solution has been found by Berman [3] for the Bertolami 
solution for the Brans-Dicke theory with a time dependent cosmological constant. Moreover, Ref.2 
predicts that G oc R for R — > 00 (asymptotically). 

8.5 An inflationary solution 

This corresponds to the case n = 1. Hence, eq.(8.17) gives 

H = C = const.. (8.33) 

Therefore 

R = const. exp(Ct) (8.34) 



84 



This solution has been obtained by [9, 10]. Hence eqs.(8.14) and (8.15) become 



p = N exp(-3jH t) , N const 



(8.35) 



and 



G = Mexp[(37 - 2)H t], M const.. 



(8.36) 



8.6 Concluding remarks 



We have analyzed a flat viscous cosmological model with varying G and A. The gravitational constant is 
shown to increase quadratically with time in the pure radiation era and linearly in the matter dominated 
era. We have also shown that Abdel Rahman model approaches this model asymptotically, i.e. for 
t 3> Rq. The model assumes to solve the horizon and monopole problems of the standard model. In this 
model we relax the assumption of the critical density. The cosmological constant retains its evolution 
with time, A oc t~ 2 , since this is assumed to be fundamental. We see that the viscosity becomes more 
important in the matter than in the radiation epoch. This model corresponds to one of the models we 
have considered recently. 
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8.8 Appendix A 



fj In order to obtain the equations describing the motion of a viscous fluid, we have to include some 



additional terms in the equation of motion of an ideal fluid. The equation of motion of a viscous fluid 
may therefore be obtained by adding to the ideal momentum flux a term a' ik which gives the irreversible 
viscous transfer of momentum in the fluid. Thus we write the momentum flux density in a viscous fluid 
in the form 



The general form of the tensor a' ik can be found as follows. Processes of internal friction occur in a 
fluid only when different fluid particles move with different velocities, so that there is a relative motion 
between various parts of the fluid. Hence a' ik must depend on the space derivatives of the velocity. If 
the velocity gradients are small we may suppose that the momentum transfer due to viscosity depends 
only on the first derivatives of the velocity. To the same approximation, a' ik may be supposed a linear 
function J^-. There can be no term in a' ik independent of J-^-, since o' ik must vanish for v = constant. 
Next, we notice that a' ik must also vanish when the whole fluid is in uniform rotation, since it is clear 
that in such a motion no internal friction occurs in the fluid. The sum 



are linear combination of the derivatives and vanish when v = !lxr, where £1 is the angular velocity. 
Hence a' ik must contain just these symmetrical combinations of the derivatives Jjjp-. The most general 
tensor of rank two satisfying the above condition is 



with coefficients £ and rj independent of the velocity. In making this statement we use the fact that the 
fluid is isotropic, as a result of which its properties must be described by scalar quantities only ( in this 
case, C; 7 ?)- The constants £ and rj are called coefficients of viscosity. Note that £ and rj are functions of 
temperature and pressure and therefore are not constants throughout the fluid. 



T" ik = pviVk - pSik + a' ik . 



dvi dv k 



dxk dxi 




8.9 Appendix B 



The apparent luminosity is related to the absolute one as by 



LR 2 





1 See e.g. L.D. Landau and E.M. Lifshitz, Fluid Mechanics (Pergamon Books Ltd. 1987) 
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where T\ is the source coordinate, and ii,io are the time of emission and reception respectively. Hence 

- f^-, (8.38) 



But 

r*° di r 1 dr 
tl W) = Jo VI - kr 2 
The scale factor can be expanded in power of (to — t) as 



= f(ri). (8.39) 



R(t) = R(t - (t - t)) = R(t )[l - (t - t)Ho - l/2(t - t) 2 q H 2 + ....] (8.40) 

where f(r±) = sin -1 n, r±, sinri, according to whether fe = 1,0, — 1. 
The red-shift z is 

z = (t - h)H + (t - ti) 2 ((7o/2 + l)H 2 + .... (8.41) 
Inverting this, the look-back time is 

to - h = H Q l z - H^(l + qo/2)z 2 + ... (8.42) 

Therefore 

n + 0(r 3 ) = i? ^to - h + l/2ff (to " *i) 2 + - 
= (iJo^o)- 1 ^- 1/2(1 +qo)z 2 + ..] 
d L = Hq^z + 1/2(1 - q )z 2 + ...] (8.43) 

Hence 

£ is usually expressed in terms of apparent bolometric magnitude m. This is defined by 

I = lO" 2 ™/ 5 x 2.52 x 10" 5 erg/cm 2 /s. 

The absolute luminosity is defined as the apparent magnitude of the source would have at a distance of 
lOpc. 

L = 1(T 2M / 5 x 3.02 x 10 35 erg s _1 . 
The distance modulus can be found as follows 

d L = itf+i™-M)/s pc 

and hence 

m - M = 25 - 51ogF + 1-086(1 - q )z + ... 

or 

m — M = 5 log cLl — 5 



